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S ■ ABSTRACT 

^ ' A new classification of neutron star cooling scenarios, involving either "minimal" cooling or 

"enhanced" cooling is proposed. The minimal cooling scenario replaces and extends the so-called 
standard cooling scenario to include neutrino emission from the Cooper pair breaking and forma- 
, tion process. This emission dominates that due to the modified Urea process for temperatures 

' close to the critical temperature for superfluid pairing. Minimal cooling is distinguished from 

enhanced cooling by the absence of neutrino emission from any direct Urea process, due either 
to nucleons or to exotica such as hyperons, Bose condensates or deconfined quarks. Within the 
minimal cooling scenario, theoretical cooling models can be considered to be a four parameter 
family involving the equation of state (including various compositional possibilities) of dense 
matter, superfluid properties of dense matter, the composition of the neutron star envelope, and 
the mass of the neutron star. The consequences of minimal cooling are explored through ex- 
tensive variations of these parameters. The results are compared with the inferred properties of 
thermally-emitting neutron stars in order to ascertain if enhanced cooling occurs in any of them. 

All stars for which thermal emissions have been clearly detected are at least marginally consis- 
tent with the lack of enhanced cooling, given the combined uncertainties in ages and temperatures 
or luminosities. The two pulsars PSR 0833-45 (Vela) and PSR 1706-44 would require enhanced 
cooling in case their ages and/or temperatures are on the lower side of their estimated values 
whereas the four stars PSR 0656+14, PSR 1055-52, Gcminga, and RX J0720.4-3125 may require 
some source of internal heating in case their age and/or luminosity are on the upper side of their 
estimated values. The new upper limits on the thermal luminosity of PSR J0205-f 6449 (in the 
supernova remnant 3C58) and RX J0007. 0-1-7302 (in CTA 1) are indicative of the occurrence of 
some enhanced neutrino emission beyond the minimal scenario. 



Subject headings: Dense matter — equation of stat^ — neutrinos — stars: neutron 



1. INTRODUCTION 

Within the last several years, several candidates 
for thermally-emitting neutron stars have been 
discovered (see, e.g., Pavlov & Zavlin (2002) for 
a short review). These stars are presumably cool- 
ing through the combination of neutrino emission 
from the interior and photon cooling from the sur- 
face, the latter responsible for their observed ther- 
mal emissions. Their temperatures have been de- 
duced by fitting atmosphere models to their spec- 
tra whereas ages can be inferred from kinematics 
or from the associated pulsar spin-down timescale. 
Neutron star cooling depends upon the equation 
of state (EOS) of dense matter as well as the neu- 
tron star mass and their envelope composition. It 
has been hoped that comparing theoretical cooling 
curves, i.e., the temperature-age or luminosity-age 
relation for neutron stars, with observations could 
yield information about their internal properties. 

A goal of this paper to explore how observations 
of thermal emission from neutron stars might be 
able to constrain the equation of state of dense 
matter. We will treat theoretical neutron star 
cooling trajectories as a four parameter series of 
models. The parameters are: 

(1) the equation of state (including various com- 
positional possibilities), 

(2) superfluid properties of the relevant compo- 
nents, 

(3) the envelope composition, and 

(4) the stellar mass. 

In nature, only one equation of state and one 
set of superfluid properties is realized, but at the 
present time, the theoretical range of superfluid 
properties for a given equation of state is so broad 
that these must be treated as an independent pa- 
rameter. In the future, as more observations be- 
come available, it should be possible to eliminate 
some combinations of EOS and superfluid param- 
eter sets, if not entire families of possibilities. It is 
also not known if a neutron star's envelope com- 
position is unique (at least for a given mass) or if 
it varies from star to star or as a fmiction of time. 
The neutron star mass is constrained to lie be- 
tween the maximum mass (a function of the equa- 
tion of state) and a minimum mass of about 1.1 
M© (set by theoretical considerations of neutron 



star birth (sec, e.g.. Burrows & Lattimer (1986)). 

Historically, theoretical neutron star cooling 
models have fallen into two categories, "standard" 
cooling or enhanced cooling. The so-called "stan- 
dard cooling" scenario has no "enhanced cooling" 
which could result from any of the direct Urea pro- 
cesses involving nucleons, hyperons, meson con- 
densates or quark matter (see, e.g. Pethick (1992) 
and Prakash (1998)). Until recently, "standard" 
cooling has been treated as being dominated by 
the modified Urea process (Friman & Maxwell 
1979). However, in the presence of superconduc- 
tivity or superfluidity in the neutron star interior, 
an additional source of neutrino emission. Cooper 
pair breaking and formation, occurs (Flowers, Ru- 
derman & Sutherland 1976; Voskresensky & Sena- 
torov 1987). For temperatures near the associated 
gap energies. Cooper pairs, in fact, dominate the 
neutrino emissivities. Although the magnitude of 
the superfluid gap energies as a function of density 
is somewhat uncertain at present, it is generally 
accepted that superfluidity occurs in neutron star 
matter. For this reason, we embark on the "mini- 
mal cooling" scenario in which "standard cooling" 
is extended to include the effects of superfluidity, 
including Cooper pair breaking and formation. 

The purpose of this paper is to explore in as 
complete a fashion as possible the consequences 
of this minimal cooling paradigm, employing the 
four kinds of parameters described above, and to 
compare our results with the inferred properties 
of cooling neutron stars. In this way, it will be- 
come apparent to what extent one or more of the 
so-called enhanced cooling mechanisms might be 
necessary to understand the observations. A fu- 
ture paper will explore in a similar fashion the 
consequences of enhanced cooling. 

One consequence of the minimal cooling paradigm, 
that enhanced cooling will not occur, is a restric- 
tion upon the equation of state involving the sym- 
metry energy. It is well-known (see, e.g., Lattimer 
et al. (1991)) that the density dependence of the 
nuclear symmetry energy controls the charge frac- 
tion in uniform beta-equilibrium matter. Since 
the direct Urea process occurs in uniform beta- 
equilibrium matter when the charge fraction ex- 
ceeds 1/9 (in the absence of muons or hyperons), 
minimal cooling thus restricts the density depen- 
dence of the nuclear symmetry energy: the crit- 
ical density for the onset of the direct Urea pro- 
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cess must remain above the star's central density. 
When muons are considered, this critical density, 
for a given equation of state, is slightly lowered. 
The appearance of hyperons can also trigger other 
direct Urea processes (Prakash et al. 1992). Thus, 
to the extent that minimal cooling can explain ex- 
isting observations, a constraint on the equation 
of state could be inferred. 

In § 2, observations of cooling neutron stars are 
reviewed. The input physics, including the equa- 
tion of state, superfluid properties, and neutrino 
emissivities are discussed in § 3. The influence of 
the neutron star envelope is briefly discussed in 
§ 4. The results of cooling calculations for min- 
imal cooling models are extensively discussed in 
§ 5. The coldest stars possible within the mini- 
mal cooling scenario are identifled in § 6 and § 7 
contains a summary of the confrontation of the 
minimal cooling paradigm with existing data. A 
comparison with other studies is performed in § 8. 
Conclusions are offered in § 9. 

2. DATA ON COOLING NEUTRON 
STARS 

Observations of neutron stars whose thermal 
emission has been unambiguously detected give 
rise to the information summarized in Tables 1 and 
2 whereas Table 3 contains results about objects 
for which only upper limits have been set. Ta- 
bles 1 and 2 display four inferred quantities: the 
total thermal luminosity L^o, the surface temper- 
ature Too, the distance d and the age t, whereas 
Table 3 omits Too • The subscript oo refers to quan- 
tities observed at the Earth which are redshifted 
relative to their values at the stellar surface. The 
data in these tables are taken from references that 
are detailed in Appendix A. In cases where a 
range of estimates is presented in these references, 
the particular parameters selected for inclusion in 
these tables are elaborated in Appendix A. Ta- 
ble 1 presents properties as inferred from models 
incorporating atmospheres dominated by hydro- 
gen, whereas Table 2 presents properties inferred 
from blackbody or heavy-element dominated at- 
mospheres. The stars displayed in Table 1 are a 
subset of those in Table 2 because the inferred 
radii of the excluded stars are far above theoret- 
ically plausible values for neutron star radii (see 
discussion below). The objects listed in Table 3 



are fainter than those listed in the first two ta- 
bles and the upper limits listed have been obtained 
only very recently thanks to the extended capabil- 
ities of XMM-Newton and Chandra. As a result, 
much less information has been obtained and the 
data analysis has not been as detailed as for the 
stars in Table 1 and 2. Moreover, in four cases no 
compact object has been detected. Indicated by a 
"?" in Table 3, these compact remnants may con- 
tain isolated black holes instead of neutron stars. 

For the purposes of restricting cooling models, 
we will not use two of these sources. RBS 1223 
is suspected of being a magnetar, judging from its 
inferred extremely high magnetic field. In addi- 
tion, its period derivative is highly uncertain. RX 
J0720. 4-3125 has a very imcertain age estimate, 
and both objects possibly have additional heat- 
ing sources compared to the other sources. RX 
J0720.4-3125 will be included in the relevant fig- 
ures, however, for purposes of comparison. 

2.1. Temperatures 

The estimation of Too and Too from the ob- 
served spectral fluxes requires atmospheric model- 
ing in which three additional factors are involved: 
the composition of the atmosphere, the column 
density of x-ray absorbing material between the 
star and the Earth, and the surface gravitational 
redshift (the surface gravity does not play a ma- 
jor role in fitting broad spectral fiux distributions. 
The column density is important because the bulk 
of the emitted flux from neutron stars is absorbed 
(mostly by interstellar hydrogen) before it reaches 
the Earth. The surface gravitational redshift, al- 
though not a factor in blackbody models, can 
influence heavy-element atmosphere models. In 
many references, the gravitational redshift was not 
optimized, but was set to the canonical value 0.304 
implied by M = 1.4 Mq and i? = 10 km. 

Since narrow spectral lines are not observed 
in any of the stars in Tables 1 and 2, the at- 
mospheric composition of these neutron stars is 
unknown. However, some information can be de- 
duced from the shape of the spectral distribution. 
Broadly speaking, neutron star atmospheres can 
be described as being either light-element (i.e., 
H or He) or heavy-element dominated. Heavy- 
element atmospheres have spectral distributions 
more closely resembling the blackbody distribu- 
tion than do light-element atmospheres (Romani 
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Table 1 

Neutron Star Properties with Hydrogen Atmospheres 



Star 


logio isd 


logio tkin 


logio Too 




logic 


Loo 


yr 


yr 


K 


kpc 


erg 


:/s 


RX J0822-4247 


3.90 


q r7+0.04 
-0.04 


0-^^-0.04 


1.9-2.5 


33.85 - 


- 34.00 


IE 1207.4-5209 


r cq+0.44 

o.oo_o 19 


n OC-H0.48 

>>-oO_0.48 


6 21+°-°'^ 

"•^-■--0.07 


1.3-3.9 


33.27- 


- 33.74 


RX J0002+6246 




o OR+O.OB 

-^-yo-o.os 


6 03+° °^ 
o-^J-o.os 


2.5-3.5 


33.08 - 


- 33.33 


PSR 0833-45 (Vela) 


4.05 


^•^O_0.31 




0.22-0.28 


32.41 - 


- 32.70 


PSR 1706-44 


4.24 




c O+0.13 
"3-O-0.13 


1.4-2.3 


31.81 - 


- 32.93 


PSR 0538+2817 


4.47 




D-UO_o.io 


1.2 


32.6- 


-33.6 



References are cited in Appendix A 



Table 2 

Neutron Star Properties with Blackbody Atmospheres 



Star 


logio isd 


logic tkin 


logio ^^oo 




d 


logio Loo 


yr 


yr 


K 


km 


kpc 


erg/s 


RX J0822-4247 


3.90 


q r-7+0.04 
-J-O' -0.04 


6.65^;?:°1 


1 - 1.6 


1.9-2.5 


33.60 - 33.90 


IE 1207.4-5209 


c cq-l-O.44 


q oc+0.48 


6.48t-;5} 


1.0-3.7 


1.3-3.9 


32.70 - 33.88 


RX ,J0002-f 6246 




"-•■^"-0.08 




2.1 - 5.3 


2.5-3.5 


32.18 - 32.81 


PSR 0833-45 (Vela) 


4.05 


4 26+° " 


o -1 n+0.02 

b.l8_o Q2 


1.7-2.5 


0.22 - 0.28 


32.04 - 32.32 


PSR 1706-44 


4.24 




99+0.04 
"■^^-0.04 


1.9-5.8 


1.8-3.2 


32.48 - 33.08 


PSR 0656+14 


5.04 




r 7-1+0.03 

^- ' -"--0.04 


7.0-8.5 


0.26-0.32 


32.18 - 32.97 


PSR 0633+1748 (Geminj 


ja) 5.53 




r 7r+0.04 


2.7-8.7 


0.123-0.216 


30.85 - 31.51 


PSR 1055-52 


5.43 




r Q2+0 02 


6.5-19.5 


0.5-1.5 


32.07-33.19 


RX J1856. 5-3754 




c yn+O.OS 


5.6-5.9 


> 16 


0.105-0.129 


31.44-31.68 


RX J0720.4-3125 


6.0 + 0.2 




5.55 - 5.95 


5.0 - 15.0 


0.1 - 0.3 


31.3-32.5 



References are cited in Appendix A 
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Table 3 

Properties of Barely Detected or Undetected Objects 



Star (SNR) ^ "^"^ 

^ yr yr kpc erg/s 



CXO J232327.8+584842 (Gas A) 




2 


.51 


3.3 


-3.6 




< 34.5 


J0205+6449 (3C 58) 


3.74 


2 


.91 


2.6 


- 3.2 




< 33.0 


PSR J1124-5916 (G292.0+1.8) 


3.45 


3.15 


- 3.30 


5 


- 6 




< 33.3 


RX J0007.0+7302 (CTA 1) 




4.0 


-4.2 


1.1 


- 1.7 




< 32.3 


? (G084.2-0.8) 




3.5 


-4.0 




4.5 


< 


30.68 - 31.45 


? (G093.3-6.9) 




3.3 


-4.0 


2.1 


- 2.9 


< 


30.65 - 31.55 


? (G127.1+0.5) 




3.3 


-3.9 


1.2 


- 1.3 


< 


29.6 - 30.75 


? (G315.4-2.3) 




3.5 - 


-4.17 


2.4 


- 3.2 


< 


30.65 - 31.80 


PSR J0154+61 


5.29 






1.7 


- 2.2 




< 32.14 



References are cited in Appendix A 



1987). This is due to the higher opacities of heavy 
elements, and seems to be the case even in the 
presence of strong magnetic fields. Since the wave- 
length range of available x-ray spectra is relatively 
small, it is possible to fit x-ray spectra with both 
kinds of atmosphere models. In general, an x-ray 
spectrum that is fit with a light-clement atmo- 
sphere will predict the star to have a lower tem- 
perature and a larger angular size than will be the 
case if a heavy-element atmosphere or blackbody 
is assumed. If the distance is known, the neutron 
star radius can be inferred (see more on this in 
§ 2.3). In some cases, fitting a star with a light- 
element atmosphere results in a predicted neutron 
star radius much larger than the canonical range 
of 10-15 km. In other cases, fitting a star with a 
heavy-element atmosphere could result in an in- 
ferred radius that is too small. 

Chang & Bildsten (2003a, b) have discussed a 
trend observed from atmospheric modeling of ther- 
mal neutron star spectra (Pavlov 2000): the in- 
ferred neutron star radii for stars younger than 
about 10^ years are consistent with canonical val- 
ues only if they are modeled with light-element at- 
mospheres (magnetized or non-magnetized). Stars 
older than about 10^ years, on the other hand, 
have inferred radii close to the canonical range 
only when modeled with heavy-element atmo- 
spheres. For this reason. Table 1 is limited to stars 
with ages less than about 10^ years, and it displays 



results inferred from modeling them with H atmo- 
spheres. Table 2, on the other hand, lists all stars, 
and displays properties deduced from blackbody 
(or heavy-element dominated) models. The tem- 
peratures and luminosities are plotted in Figure 1, 
and are also selected according to this trend and 
our desire that the inferred stellar radius lies in 
a theoretically plausible range. The temperature 
and luminosity are taken from Table 2 unless val- 
ues for them appear in Table 1. 

The above trend implies that the atmospheric 
composition of a neutron star evolves from light 
to heavy elements with a timescale of about 10^ 
years. This possible evolution is considered in 
more detail in § 5.2. 

2.2. Ages 

The precise ages of observed cooling neutron 
stars are not always known. Most stars listed in 
Tables 1 through 3 are known radio and/or x- 
ray pulsars and their ages can be estimated from 
the oberved spin-down rate using tsd = P/'^P, 
where P and P are the period and its time deriva- 
tive, respectively. In some cases, kinematic infor- 
mation is available and ages can be inferred by re- 
lating pulsar transverse velocities to the distances 
from the presumed sites of origin as, e.g., the ge- 
ometric center of the associated supernova rem- 
nant or a nearby cluster of massive OB stars. In 
the case of an association with a supernova rem- 
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Log t [yrs] 

Fig. 1. — Inferred temperature Too (top) and lu- 
minosity Loo (bottom) versus age for neutron stars 
with thermal emission. Data from Table 1 are 
marked as "H atmosphere fits" and data from Ta- 
ble 2 as "Blackbody fits" . 
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PSR J0154+61 



Fig. 2. — Inferred upper limits on thermal lumi- 
nosity Loo, versus age, for the compact objects 
listed in Table 3. 



nant, the age can also be estimated by the general 
properties of the remnant and, in the best cases, 
by association with historical supernovae. We will 
generically refer to these various alternatives to 
tsd as the "kinematic age" tkin- Both ages, where 
available, are shown in Tables 1, 2 and 3. The 
data in these tables show that, in general, there is 
a large discrepancy between the spin-down age and 
the kinematic age in cases where both are given. 
In most cases, the spin-down age is longer, but in 
the case of the Vela pulsar, it is shorter. Typical 
discrepancies are of order 3 or larger. For this rea- 
son, we have used the kinematic age in Figures 1 
and 2 where available, and otherwise have assigned 
an uncertainty of a factor of 3 in each direction to 
the spin-down ages. 

2.3. Distances and Luminosities 

The distances are estimated from pulsar disper- 
sion measures, estimated distances to the related 
supernova remnants, or observations of interstel- 
lar absorption to other stars in proximity. In three 
cases, parallax estimates are available. The details 
are discussed in Appendix A. Uncertainties in the 
distances are in many cases rather large. Since the 
inferred luminosities of the stars are proportional 
to the square of the assumed distances, it is usu- 
ally the case that the inferred stellar luminosity 
has greater relative error bars than the inferred 
stellar temperature. However, in cases in which 
the composition of the stellar atmosphere is uncer- 
tain, but the distance to the source is accurately 
known, the inferred stellar luminosity might be 
more accurately estimated. A consistency check 
of the measurements of Too, L^o and the distance 
d is that the relation 



(1) 



should give a radius at infinity i?oo comparable to 
the radius of a neutron star. This is the case for 
the measurements listcid in Table 1, whereas for 
the measurements based on BE spectral fits listed 
in Table 2 only 1055-52 has a possibly acceptable 
i?oo, but with very large errors due to the uncer- 
tainty in d. However, BB models are overly sim- 
plistic. Non-magnetic heavy-element dominated 
atmospheres tend to have values of i?oo factors of 
2 to 3 larger than a BB (Romani 1987), so that 
essentially all the sources listed in Table 2 satisfy 
this consistency check. For the objects listed in 
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Table 3, this consistency test is only marginally 
possible for J0205+6449 (3C58). 

Theoretical cooling calculations also give an 
effective temperature and a luminosity L°° 
which are related to each other by the equa- 
tion (see also equation (B9)) 

L°° = e2*(^)L(i?) = 4TrR°°^ -assT^^, (2) 

where we have used the superscript oo to denote 
the theoretical values and subscript cxd for the ob- 
served values at infinity in order to emphasize the 
difference. It is only for a star for which the mea- 
sured Too and Lqo satisfy equation (2), and for 
which an accurate measurement of d exists (im- 
plying a small error bar on Loo), that comparison 
of cooling curves with data in terms of T or L are 
equivalent. For the stars listed in Table 2 which 
do not pass the above consistency test, the mea- 
sured Too is thus not an effective temperature and 
cannot be directly compared with the calculated 
Tg°°. In these cases, the luminosity Too is more 
representative of thermal emission and should be 
used for comparison with L°°. For this reason, 
we have chosen to tabulate luminosities as well as 
temperatures in Tables 1 and 2, and have plot- 
ted both temperature and luminosity in Figure 1. 
In Table 3, wc have reported only upper limits to 
Too, which is the quantity that observation can 
usefully constrain, and plotted them separately in 
Figure 2. 

One feature notable in Figures 1 and 2 is the 
sizes of the error boxes, particularly in the age di- 
mension. These uncertainties represent an inher- 
ent difficulty in using these observations to firmly 
constrain the details of neutron star cooling. For 
this reason, instead of attempting to detail prop- 
erties of the equation of state, superconductiv- 
ity, and/or neutrino emissivities from the observa- 
tions, our approach will be to model a reasonably 
broad range of acceptable physical inputs in order 
to determine ranges of parameters that might be 
excluded by the present data. 

3. INPUT PHYSICS 

The standard general relativistic equations de- 
termining the structure and thermal evolution of 
a neutron star are briefly summarized in Ap- 
pendix B. Given an equation of state (EOS), de- 
scribed in § 3.1, we solve numerically the TOV 



equations of hydrostatic equilibrium, and build 
our stars. The equations of energy conserva- 
tion, equation (B2), and energy transport, equa- 
tion (B5), with their corresponding boundary con- 
ditions are then solved numerically with a fully 
general relativistic Henyey-type stellar evolution 
code specially developed for neutron stars (Page 
1989). The required physics input are described 
in the next section. The outer boundary condi- 
tion, equation (B6), is implemented in terms of an 
envelope, described in § 4. 

3.1. The Equation of State 

The gross properties of a neutron star (such 
as its mass and radius) and its interior compo- 
sition (which influences the thermal evolution) 
chiefly depend on the nature of strong interac- 
tions in dense matter. Investigations of dense mat- 
ter can be conveniently grouped into three broad 
categories: nonrclativistic potential models, effec- 
tive field theoretical (EFT) models, and relativis- 
tic Dirac-Brueckner-Hartree-Fock (DBHF) mod- 
els. In addition to nucleons, the presence of soft- 
ening components such as hyperons, Bose con- 
densates or quark matter, can be incorporated in 
each of these approaches. Some general attributes, 
including references and typical compositions, of 
equations of state (EOS's) in each of these ap- 
proaches have recently been summarized by Lat- 
timer & Prakash (2001). 

In this work, we employ four EOS's in the cat- 
egory of nonrclativistic potential models in which 
only nucleonic degrees of freedom are considered. 
Two of these are taken from the calculations of 
the Argonne and Urbana groups. The EOS la- 
beled WFF3, from Wiringa, Fiks & Fabrocini 
(1988), is based on the variational calculations us- 
ing UV14+TNI potential and that labeled APR, 
from Akmal, Pandharipande & Ravenhall (1998), 
utilizes the AV18 potential plus the UIX potential 
plus the Svf, boost. APR represents the most com- 
plete study to date of Akmal & Pandharipande 
(1997), in which many-body and special relativis- 
tic corrections are progressively incorporated into 
prior models including that of WFF3. 

For isospin symmetric matter, the equilibrium 
densities of the WFF3 and APR models are no = 
0.163 fm^'^ and 0.16 fm~^, respectively, with cor- 
responding compression modulii of 269 MeV and 
274 MeV, respectively. In isospin asymmetric 
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matter, the density dependent symmetry energy 
S{nb,x) is defined by the relation 

E{nb,x) = E{nb,l/2) + S{nb,x), (3) 

where E is the energy per particle, nb = Un + rip 
is the baryon number density, and x = Up/ub is 
the proton fraction. In practice, S{nb,x) can be 
expanded as 

S{nb, x) = S2{nb){l - 2xf + 54(n,,)(l - 2x)^ ■■■ , (4) 

where the term involving ^4 is generally very 
small. S{nb,x) plays a crucial role in a neutron 
star's thermal evolution insofar as it determines 
the equilibrium proton fraction, which in turn de- 
termines whether or not the direct Urea process, 
n — > p + e~ + Pe, is permitted to occur in charge 
neutral beta-equilibrated matter (Lattimer et al. 
1991). 

The equilibrium proton fraction is determined 
from the condition 

He = /i = - = -{dE/dx) (5) 
= 4(1 - 2x)[S2{nb) + 2S4{nb){l - 2xf + •••], 

where /i^ {i = e,n,p) are the chemical poten- 
tials. For ultrarelativistic and degenerate elec- 
trons, He = tlc{3n'^nbxy^^, since due to charge 
neutrality rie = Up in matter in which the only 
leptons are electrons. 

When the electron Fermi energy is large enough 
(i.e., greater than the muon mass), it is ener- 
getically favorable for the electrons to convert to 
muons through e~ IJ,~ -\- + Ve ■ Denoting 
the muon chemical potential by /i^, the chemical 
equilibrium established by the above process and 
its inverse is given by ^ti = Me- At the thresh- 
old for muons to appear, /i^ = ra^c? = 105 MeV. 
Noting that the proton fraction at nuclear den- 
sity is small, one has the approximate relationship 
AS2{u) /ra^c^ ^ 1, where u — fib/n^). Using a 
typical value S2{u = 1) ~ 30 MeV, one may ex- 
pect muons to appear roughly at nuclear density 
no = 0.16 fm~^. Above the threshold density, 

/Ue = A = A*M = \J {hcyiSn'^nbX^)'^/^ + ml^d^ , (6) 

where = Ufj^/ub is the muon fraction in mat- 
ter. The charge neutrality condition now takes 
the form ne+rifj^ = rip, which, together with the 



relation of chemical equilibrium in equation (6), 
establishes the lepton and proton fractions in mat- 
ter. The appearance of muons has the consequence 
that the electron fraction .t,. = n^/nb is lower than 
its value without the presence of muons. 

If S{nb, x) does not rise sufficiently rapidly with 
density, the equilibrium proton fraction will re- 
main below a critical value (of order 11% in mat- 
ter with e~ only and 14% in matter with both 
e~ and /i~) required for the direct Urea process. 
The critical proton fraction is determined by re- 
quiring simultaneous conservation of energy and 
momentum among the participating fermions. In 
this case, cooling occurs via the modified Urea pro- 
cess, n + n ^ n + p + e~ + De, and several other 
similar processes (see § 3.5), modulated by effects 
of possible nucleon superfluidity. Because of the 
additional fermions involved, the emissivity of the 
modified Urea process is several orders of magni- 
tude lower than the direct Urea process. 

The symmetry energies 52(no) of the WFF3 
and APR models at their respective equilibrium 
densities no have the values 29.5 MeV and 32.6 
MeV, respectively. In the WFF3 model, the sym- 
metry energy 52 (nb) rises slowly with density and 
X never reaches the critical value for the direct 
Urea process. However, in the case of the APR 
model, the direct Urea process becomes possible 
at nB > 0.78 fm~'^, which corresponds to a neu- 
tron star mass of Mcr = 1.97 Mq . For this model, 
we will therefore consider only stars with masses 
below this threshold. 

We also consider two EOS's from the phe- 
nomenological non-relativistic potential model of 
Prakash et al. (1997) which is designed to repro- 
duce the results of more microscopic calculations 
at zero temperature, and which allows extensions 
to finite temperature. The EOS's chosen are la- 
beled BPAL21 and BPAL31, which have bulk nu- 
clear matter incompressibilities Kg = 180 or 240 
MeV, respectively. In both eases, the symmetry 
energy, at the empirical symmetric matter equi- 
librium density of no = 0.16 fm~^, was chosen 
to be 30 MeV. Furthermore, the potential part of 
S{nb,x) varies approximately as ^/ub/no in both 
cases, which is close to the behavior exhibited in 
the EOS of APR. 

In Figure 3, the symmetry energies (top 
panel) and corresponding proton fractions (bot- 
tom panel) in charge-neutral beta-stable neutron 
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star matter are shown. In Figure 4, the pressure 
of neutron star matter is shown as a function of 
baryon density for the EOS's considered in this 
work. The differences in the high-density behav- 
ior of these two EOS's are largely attributed to 
differences in the underlying three-body interac- 
tions. 

The reason why we do not consider EOS's based 
on effective field-theoretical (EFT) and relativis- 
tic Dirac-Brueckner-Hartree-Fock (DBHF) models 
in this work merits some discussion. In EFT ap- 
proaches based on the prototype Walecka model, 
interactions between nucleons are mediated by the 
exchange of cr— , ui—, and p— mesons. At the mean 
field level, the symmetry energy in this approach 
is given by (Horowitz & Piekarewicz (2001)) 

^^^"^^ = 6Vfc^ + M-+8(i| + 2/(.o,a;o))' ' 

where ctq, lvq, and po are the mean- field expec- 
tation values of the fields, g^- and Qp are the 
a— and p— meson couplings to the nucleon, and 
M* = M — gcrCTo is the nucleon's Dirac effective 
mass. The quantity /(cro,wo) summarizes effects 
of density dependent nonlinear interactions aris- 
ing from a—, oj—, and p— mixings and have re- 
cently been employed to explore deviations from 
the linear behavior of the second term with den- 
sity in the case / = 0. When the symmetry en- 
ergy rises linearly with density, the critical proton 
fraction for the direct Urea process is reached at 
2 — 3 no, which is well within the central densities 
of both 1.4M0 and maximum mass stars obtained 
with these EOS's. It is possible, however, to forbid 
the direct Urea process with a suitable choice of 
/ ^ (Steiner et al. 2004). These eases, however, 
resemble the potential models considered above. 

All DBHF calculations reported thus far in the 
literature (e.g. Miither, Prakash & Ainsworth 
(1987); Engvik et al. (1994)) find that proton frac- 
tions favorable for the direct Urea process to occur 
are reached in stars whose masses are larger than 
~ 1.3Mq. Since our intention here is to explore 
the extent to which model predictions can account 
for observations without invoking the direct Urea 
process and its variants involving hyperons, Bosc 
condensation or quarks, we defer a discussion of 
these models to a separate work. 

In Figure 5, we show the mass versus radius and 



versus central density curves for the four EOS's 
chosen. Features of relevance to the discussion of 
cooling to note are: 

• The radii of maximum mass configurations 
(1.7 < Mmaxf^Q < 2.2) are confined to the 
narrow range 9-10 km. 

• The radii of 1.4M0 stars lie in the narrow 
range 11 — 12 km. 

Significant deviations from such a tight clustering 
of radii occur only in those eases in which 

• normal nucleonic matter is described through 
the use of EFT or DBHF models. Lattimer 
& Prakash (2001) showed that the neutron 
star radius is proportional to the density 
derivative of the symmetry energy in the 
vicinity of nuclear matter density. There- 
fore, in this case, relatively large radii for 
both 1.4M0 and maximum mass configu- 
rations occur. Similarly, a relatively large 
density dependence of the symmetry energy 
also permits the direct Urea process to occur 
in this case. 

• extreme softening is induced by the pres- 
ence of additional components such as hy- 
perons, Bose condensates or quarks: In this 
case, significantly smaller radii are possi- 
ble. Such components also lead to relatively 
rapid cooling. 

As mentioned earlier, both of these cases fall out- 
side the "Minimal Cooling Scenario" and will be 
investigated separately. 

For completeness, we note that in all four eases 
considered, the crust-core transition occurs at 
nb ~ 0.1 fm~^ or equivalently p ~ 1.6 x 10^^ 
g em~^. For the EOS in the crust region, we 
employ the EOS of Negele & Vautherin (1973) 
above neutron drip and that of Haensel, Zdunik, 
& Dobaczewski (1989) below neutron drip. 

3.2. Nucleon Effective Masses 

Under degenerate conditions (T ^ p) and in 
the absence of collective excitations close to the 
Fermi surface, physical quantities such as the spe- 
cific heat, entropy, and superfluid gaps, and pro- 
cesses such as neutrino emission from particles 
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Fig. 3. — Symmetry energy (top panel) and proton 
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Fig. 4. — Pressure vs baryon density for the four 
EOS's employed in this work. 
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with energies in the neighborhood of the Fermi 
energy, depend sensitively on the so-called Landau 
effective masses of particles. Formally, the Landau 
effective mass m* of any degenerate fermion (n, p, 
e, and fj. in our case) is defined by 



PF 



de{p) 




dp 


P=PF_ 



(8) 



where e{p) is the single-particle energy of the par- 
ticle with momentum p, and the derivative is eval- 
uated at the Fermi momentum pp. For nonrel- 
ativistic interacting nucleons, the single particle 
energies of the neutron and proton can be written 
as 



filled circles in Figure 7 show the symmetric mat- 
ter Landau effective masses for this case. Lack- 
ing further input, and encouraged by the fact that 
the APR results for neutron and proton effec- 
tive masses in beta-stable isospin asymmetric mat- 
ter bracket the WFF results in isospin symmetric 
matter, we take the results of the APR model to 
apply for the WFF model as well. As we will show, 
our final results are not significantly affected by 
this approximation. 

The neutron and proton Landau effective 
masses for the BPAL21 and BPAL31 models can 
be obtained straightforwardly from the single- 
particle potentials given in Prakash et al. (1997). 
Explicitly, 



P 

en{p) = 2^ + ^"('^«',P) 

ep{p) = — + Up{nt,p) (9) 

where m is the nucleon mass in vacuum, and Un 
and Up are the neutron and proton single-particle 
momentum-dependent potentials which are ob- 
tained by appropriate functional differentiations 
of the potential energy density. 

From the Hamiltonian density in Appendix A 
of APR, the neutron and proton Landau effective 
masses are 



m 



j=l,2 



{2zu 



\2/3' 



_ 2 {2Q + AZi) 
~ 5 ElR^ 



2 (C, - SZi) 



E%Rl 



(12) 



where is the Fermi energy at the equilibrium 
density no and the parameters Ci.,Zi and Ri can 
be found in Tables 1 and 2 of Prakash et al. (1997). 
The dash dotted curves in Figure 7 show the ef- 
fective masses for the BPAL models. 

The general trends to note in Figure 7 are 



m 
m 



2m 

l + ^{ps + zp5)nte-P^^^ 



(10) 



where z = (1 — x) for neutrons and z ^ x for pro- 
tons, and P3 = 89.8 MeV fm^, = 0.457 fm^, 
and p5 = —59 MeV fm^. The solid curves in Fig- 
ure 7 show the variation of the effective masses 
with density for both neutrons and protons in 
charge-neutral beta-stable matter corresponding 
to the EOS of APR. 

To date, single-particle energies for the WFF3 
model are available only for symmetric nuclear 
matter up to rib = 0.5 fm~^ (Wiringa (1988)). 
Using the parameterization from equation (7) of 
Wiringa (1988), we find that 



m 



tfiA-^^\ A2 



(11) 



where the density-dependent parameters (3 and A 
are tabulated in Table I of Wiringa (1988). The 



• the steady decrease of m* and m* with den- 
sity, and 

• m* > m* in charge- neutral beta-stable mat- 
ter. 

It is interesting to observe that there is more 
spread in the model predictions for m*, partic- 
ularly with increasing density, than for m*. 

For noninteracting relativistic particles, such as 
e and fi, the Landau effective masses are given by 



(13) 



where m denotes the appropriate vacuum mass. 

3.3. Pairing 

The Fermi surface of a degenerate system of 
fermions, as are e's, /u's, p's and n's in a catalyzed 
neutron star, becomes unstable in the presence 
of an attractive interaction between the particles 
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Fig. 6. — Fermi momenta of neutrons (upper 
curves) and protons (lower curves) versus baryon 
density for the APR and BPAL EOS's (BPAL21 
and BPAL31 give nearly the same values). 
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Fig. 7. — Ratio of the Landau effective mass m* 
to the vacuum mass m for neutrons and protons 
vs. density for the four EOS's. 



whose momenta lie close to the Fermi momentum: 
this is the Cooper theorem (Cooper 1956). As a 
result of this instability, the ground state of the 
system is reorganized with a gap A in the energy 
spectrum around the value of the Fermi energy; 
no particle can have an energy between Ep — A 
and Ef + A. 

This instability usually appears as a second or- 
der phase transition, the gap acting as the or- 
der parameter of the transition which has a corre- 
sponding critical temperature Tc. For T > Tc, the 
system behaves as a normal Fermi liquid whereas 
for T i^Tc the gap grows in magnitude which re- 
sults in superfluidity or superconductivity. The 
precise value of Tg depends on the nature of the 
pairing interaction. 

There is no obvious attractive interaction be- 
tween e's and/or /x's in a neutron star and thus 
they are not expected to become paired at tem- 
peratures relevant for our concerns here (Baym & 
Pethick 1975). 

For nucleons, the strong interaction provides 
several channels in which pairing is possible with 
Tc's of order MeV. Nucleon-nucleon scattering 
data in vacuum indicate that at low momentum 
pairing should occur through Cooper pairs with 
zero angular momentum i in a spin-singlet state, 
^5o, whereas for larger momenta an i = 1 spin- 
triplet {J = L + S = 2),'^P2 pairing becomes fa- 
vorable. Starting from a knowledge of the nucleon- 
nucleon interaction in vacuum, the difficulty of 
obtaining reliable values of the gap in a medium 
is illustrated by considering the result for the so- 
lution of the gap equation in the so called BCS 
weak coupling approximation (Bardeen, Cooper & 
Schrieffer 1957): 

^ Ef e-V^^(o) , (14) 

where V is the in-medium pairing interaction in 
the corresponding channel and A^(0) = m*pF/T^^h^ 
is the density of states at the Fermi surface. Small 
variations of to* and medium effects on V affect 
A in an exponential manner. 

The best-studied case is the n ^Sq gap in pure 
neutron matter. The gap appears at densities of 
order no or lower, where to* is well determined 
and the pairing interaction in vacuum for the 
channel is accurately known. Simple arguments 
(Pines 1971) show that medium polarization, the 
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dominant medium effect on V, should induce a re- 
duction of the ^5*0 gap (Clark et al. 1976) from its 
value without medium polarization. Much effort 
has been dedicated to take into account medium 
polarization at various levels of approximation. 
With time and improving many-body techniques, 
the results arc beginning to show a convergence 
for the maximum value of T^, which is in the range 
~ 0.5 to 0.7 X 10^° K, as can be seen in Figure 8. 
The density range in which this gap is non-zero is 
still somewhat uncertain and corresponds to the 
inner part of the crust and, possibly, the outer- 
most part of the core. 

Since the results shown in Figure 8 are for uni- 
form pure neutron matter, they will be altered by 
the presence of a small fraction of protons in the 
outer core and the nonuniformity of neutron den- 
sity due to nuclei (or nuclear clusters) in the in- 
ner crust. This latter effect has been studied re- 
cently by Barranco et al. (1997) who show that it 
does not alter significantly the results, at least at 
the level of accuracy required for the study in the 
present paper. 

The p ^Sq gap is similar to the n ^Sq gap and 
occurs at similar Fermi momenta kp, but since 
protons represent only a small fraction of the nu- 
cleons, this translates to high densities which al- 
lows the gap to persist in much deeper regions of 
the core than the n ^Sq gap. The values of Tc 
from several calculations are shown in Figure 9. 
An essential immediate difference compared to the 
n ^Sq gap is that the p ^Sq gap is much smaller, 
m* being smaller for protons than for neutrons 
(see Figure 7). It should be noted that all cal- 
culations shown in this figure have employed val- 
ues of m* larger than the values we report in Fig- 
ure 7. Insofar as the results of Figure 7 for APR 
are indicative of the likely magnitudes of m*, the 
values of Tc in Figure 9 are likely overestimated, 
particularly at large kp- Moreover, medium po- 
larization effects on V are much more difficult to 
take into account for the p ^Sq gap than for the 
neutron gap. Such effects are expected to reduce 
the size of the gap and, to date, only two works 
have attempted to include them (Niskanen & Sauls 
1981; Ainsworth, Wambach & Pines 1991). The 
estimates of Ainsworth, Wambach & Pines (1991) 
show that medium polarization reduces the 
gap roughly by a factor of three in the stellar core. 
It is important to notice that all these calculations 




kp(n) [fm-^] 

Fig. 8. — Neutron ^Sq pairing critical temper- 
ature Tc vs neutron fermi momentum kp from 
the calculations of Ainsworth, Wambach & Pines 
(1989) (labeled as "AWP 11" & "AWP III": two 
slightly different results), Wambach, Ainsworth 
& Pines (1993) ("WAP"), Chen, et al. (1993) 
("CCDK"), Shulze et al. (1996) ("SCLBL"), and 
Schwcnk, Friman & Brown (2003) ("SFB"). The 
dotted curve shows the results of Shulze et al. 
(1996) in the case where medium polarization is 
not included. Medium polarization effects reduce 
the ^5*0 gap by about a factor three. The vertical 
dotted line shows the location of the crust-core 
boundary. 
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find that Tc vanishes for kp > 1-5 fm ^ and in 
most cases for fc^ > 1 fm~^. 

The so-called n gap actually occurs in the 
^P2 — channel, since the tensor interaction cou- 
ples channels with 6L = 2. This coupling with the 
channel increases the gap (Takatsuka 1972a). 
We present an illustrative sample of published Tc 
curves in Figure 10. The large differences among 
these curves points to the inherent difficulty in pin- 
ning down the magnitude of this gap. The gap 
possibly extends to high densities where m* is un- 
certain. In addition, the presence of a small frac- 
tion of protons has generally been ignored except 
in the work of Elgar0y et al. (1996), who found a 
reduction of the gap by a factor ~ 3 when consid- 
ering neutron-proton matter in /3-equilibrium. 

A fundamental problem, emphasized recently 
by Baldo et al. (1998), is that even the best 
modern models of the nucleon-nucleon interaction 
{in vacuum) fail to reproduce the experimental 
phase shift in the channel at laboratory ener- 
gies above 300 MeV (corresponding to the pion- 
production threshold). Translating this energy 
into an equivalent density implies that the bare 
pairing interaction is not understood at densities 
■~ I.Iuq. Moreover, medium polarization effects 
have not been included in any of the calculations 
displayed in Figure 10. Estimates of such effects 
had shown that they should strongly enhance the 
^P2 - ^F2 gap (Jackson et al. 1982), but the re- 
cent results of Schwenk & Friman (2004) indi- 
cate that the medium-induced spin-orbit interac- 
tion strongly suppresses this gap. For this reason 
we will also consider the possibility that the n^P2 
gap is vanishingly small. 

The J = 2 with L = 1 or 3 angular momen- 
tum and the tensor coupling make the gap equa- 
tion a system of coupled integral equations, one 
for each value of the magnetic quantum number 
mj and the gap is not isotropic: A = A{9, <j); T), 
where and (p are the polar angles of the mo- 
mentum k. Until recently most works had looked 
for single component solutions with \rnj\ = 0, 1 
or 2, but Zverev, Clark & Khodel (2003) have 
shown that when considering multicomponent so- 
lutions there are at least 13 distinct phases: for 7 
of them A(0, 0; T) vanishes on the Fermi surface at 
some values of {0, </>) whereas the other 6 are node- 
less. Nodeless gaps are energetically favored over 
gaps with nodes, but by a small amount [see, e.g., 
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Fig. 9. — Proton pairing critical temperature 

Tc vs proton fermi momentum kp from the cal- 
culations of Takatsuka (1973) (labeled as "T"), 
Amundsen & 0stgaard (1985a) ("AO"), Chao, 
Clark & Yang (1972) ("CCY_ms" & "CCY_ps": 
two slightly different results), Chen, et al. (1993) 
("CCDK"), Baldo et al. (1992) ("BCLL"), El- 
gar0y et al. (1996) ("EEHO"), Niskanen & Sauls 
(1981) ("NS"). Only the calculation of Niskanen 
& Sauls (1981) include medium polarization. The 
shaded region shows the estimates of Ainsworth, 
Wambach & Pines (1991) of the range of values in 
which Tc should lie due to medium polarization. 
The vertical dotted line shows the location of the 
crust-core boundary. 
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Amundsen & 0stgaard (1985b)], and may become 
disfavored in the case of very fast rotation of the 
star or in the presence of an ultra-strong magnetic 
field (Muzikar, Sauls, & Serene 1980). We will, 
in this work, assume that the ^P2 — ^F2 gaps are 
nodeless, since this maximizes the effect of pairing 
on cooling and seems energetically favored, and 
more specifically assume its angular dependence 
to be A(6', <?!); T) oc (1 + 3 cos^ 6*)^/^, corresponding 
to the pure mj = phase. 

For both the isotropic ^So and the pure mj = 
phase of the ^Pj— "^^2 gap, the critical temperature 
Tc and the T = gap are related by 

ksTc « 0.57Ao (15) 

where Aq is A(T = 0) for ^Sq and the angle av- 
eraged value of A{9,(p;T = 0)^ over the Fermi 
sphere for -^i^2 (Baldo ct al. 1992). The tem- 
perature dependence of A for these two cases has 
been calculated and fitted by simple analytical ex- 
pressions by Levenfish & Yakovlev (1994a). 

Once the energy gap A is given, the quasi- 
particle energy spectrum near the vicinity of the 
Fermi surface can be expressed as 

e(k) =Ef- \/A(k)2 + 772 foj. k<kF 

= Ef + v/A(k)2 +7^2 fo^ k>kF (16) 
where the quantity 77 is given by 

^ ~ 2m* ~ 2m* 

= (^) ■ - kF) = VF-{k- kp) (17) 
\m* J 

The effect of this change in the quasi-particle spec- 
trum is discussed below in § 3.4 and § 3.5. 

Finally, it is instructive to consider Tc for the p 
^S[) gap and the n ^'^^2 gap in terms of density 
and also in terms of the volume in the stellar core, 
once an EOS and a stellar mass has been chosen. 
(This is presented for a few cases in Figures 20 and 
21 in a 1.4 M© star built with the APR EOS.) The 
p ^S'o gap vanishes or is very small in the inner core 
for most calculations whereas the n — "^^2 gap 
is more likely to reach the stellar center with high 
values. 

3.4. The Specific Heat 

The total specific heat (per unit volume) at con- 
stant volume, Cy, receives contributions from all of 




Fig. 10. — Neutron "^^2 — pairing critical tem- 
perature Tc vs neutron fermi momentum hp- The 
three thin continuous lines show the results of: 
(1) Amundsen & Ostgaard (1985b), (2) Takat- 
suka (1972b), and (3) Elgar0y et al. (1996). The 
three dotted lines show results obtained assum- 
ing m* = m„ to illustrate the strong reduction 
due to small m*„. (0) Hoffberg et al. (1970), (1) 
Amundsen & Ostgaard (1985b), and (2) Takat- 
suka (1972b). The three thick continuous lines, 
(a), (b), and (c), bracket the results of Baldo et 
al. (1998). The vertical dotted line shows the lo- 
cation of the crust-core boundary. 
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the constituents inside the star. In the homoge- 
neous phase above nuclear density, 



E 



(18) 



-e ,jU ,n,p 



For a given species of unpaired spin ^ fermions, Cy 
is given by 



1 de 



J (27^ dT J 

2 J d^k W 



(27r)3 dfi 

(19) 



where e{k) is the single-particle spectrum, / is the 
Fermi-Dirac distribution function, and n and T are 
the chemical potential and temperature, respec- 
tively. Under the degenerate conditions of inter- 
est here, T <C /i, the contribution from the second 
term above can be safely neglected with the result 



pIf 



(20) 



In writing the rightmost relation above, the den- 
sity of states at the Fermi surface A''j(0) has been 
expressed in terms of the Landau effective mass 
m* (see equation (8)) through the relations 



Ni{0) 



kpVF ' 



de 
dk 



(21) 



where vp denotes the velocity at the Fermi surface. 

The effect of pairing interactions on the specific 
heat depends on the disposition of T with respect 
to Tg. When T reaches T^, there is a sharp in- 
crease in the specific heat due to the large fluctu- 
ations characteristic of a second order phase tran- 
sition. Subsequently, when T <C Tc, a Boltzmann- 
like suppression occurs due to the presence of a 
gap in the energy spectrum, equation (16). In 
practice, these effects are taken into account by 
using control functions that would multiply the 
unpaired values of c^; these functions have been 
calculated for nucleon pairing in both the ^Sq and 
channels by Levenfish & Yakovlev (1994a) and 
are displayed in the upper panel of Figure 13. 

The cumulative contributions to Cy are pre- 
sented in Figure 11. Once a nucleonic component 



becomes paired, its contribution will ultimately be 
suppressed when T <C Tc; but in all cases the lep- 
ton contribution will always remain. For different 
temperatures the various contributions all scale as 
T in the absence of pairing and their relative im- 
portance hence remains the same as illustrated in 
this figure. 

The suppression by pairing is illustrated in Fig- 
ure 12 for the case of neutrons at various temper- 
atures beginning with 10^ K as in Figure 11. The 
results in this figure, in agreement with those in 
Figure 13, show that beyond the initial strong en- 
hancement at T ~ Tc, c„ will be almost completely 
suppressed only when T < 0.1 to 0.2 x Tc every- 
where within the core. 

In the crustal region, contributions to the spe- 
cific heat arise from the neutron gas in the inner 
crust, the degenerate electron gas and the nuclear 
lattice. The contribution of the neutron gas is 
strongly suppressed by the n ^Sq gap- The spe- 
cific heat of the crust is smaller than that of the 
leptons in the core and is hence not very impor- 
tant, but included in all our calculations. 

3.5. Neutrino Emissivities 

Until the time that photon emission takes over 

the cooling of the star, the thermal energy of the 
star is lost from both the crustal layers and the 
core of the star chiefly by neutrino emission. The 
various neutrino emission processes that are in- 
cluded in our study are summarized below. 

In the crust of the star, we include neutrino pair 
emission from plasmons according to Haft, Raffelt 
& Weiss (1994), or, from the practically equivalent 
results of Itoh et al. (1996), and from electron- 
ion bremsstrahlung according to Kaminker et al. 
(1999). These two processes are not affected by 
nucleon pairing. We also consider neutrino pair 
emission from neutron-neutron bremsstrahlung 
and its suppression by neutron pairing. The treat- 
ment of pairing suppressions is outlined below in 
the context of similar suppressions in the core. In 
addition, we consider neutrino pair emissions from 
the formation and breaking of n Cooper pairs, 
also described further below. In the presence of 
a magnetic field, synchrotron neutrino pair emis- 
sion from electrons also occurs (Bezchastnov et 
al. 1997), but contributions from this process are 
negligible. Photo-neutrino emission and neutrino 
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pairs from e^e~ pair annihilation arc effective 
only at low density and high temperature and are 
not relevant here. 

In the core of the star, we include (1) the 
modified Urea processes and the similar nucleon 
bremsstrahlung processes with their correspond- 
ing suppressions by nucleon pairing, and (2) neu- 
trino pair emission from the formation and break- 
ing of Cooper pairs. The emissivity from the neu- 
tron branch of the modified Urea process 

n + n' — > p + n' + 1 + VI 
p + n' + 1 — > n + n' + 1/1 , 



(22) 



where I is either an electron or a muon and f ; or 9i 
is the associated neutrino or antineutrino, is taken 
from Priman & Maxwell (1979) and Yakovlev & 
Levenfish (1995). Exphcitly, 
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kpo J \ kpo 



109 K 



(23) 



where kpo = 1-68 fm~^ is a fiducial normalization 
factor. The emissivity from the proton branch of 
the modified Urea process 

n + p' — » p + p' + 1 + 
p + p' + 1 — »• n+p' + 1/1 , 



(24) 



is taken from Yakovlev & Levenfish (1995) in the 
form 



\kF0j\ ^kFpj \kFoJ\ 



rrin J\mp 
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T 



Fp 



109 K 



(25) 



In equations (23) and (25), the coefficients an, ctp, 
Pn, and Pp are of order unity and describe correc- 
tions due to the momentum transfer dependence 
of the matrix element in the Born approximation 
(a„^p) and due to non-Born corrections and strong 
interaction effects beyond the one pion exchange 
plus Landau coeSicients {Pn,p) (Friman & Maxwell 
1979; Yakovlev & Levenfish' 1995). To be specific, 
following Yakovlev & Levenfish (1995) we use 



an = ap = 1.76-0.63 (1^)' 
Pn = f3p = 0.68. 



(26) 



In addition to the above two charged current 
modified Urea processes, three neutral current 
bremsstrahlung processes 



n + n' — > n + n' + 1^1 + Pi 
n + p' — »• n + p' + i'i +iyi 
p + p' — »■ P + p' + i'i + i'i , 



(27) 



where the pairs i/iDi can be an e, /i or t neu- 
trino pair, also contribute. Their emissivities are 
(Priman & Maxwell 1979; Yakovlev & Levenfish 
1995) 
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(28) 



and 



Bremnp ^3 ^ 1.5 X 10^° ' ^ ""^ ' ^ 
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(29) 
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^Brem PP = 3 x 7.4 X 10 



kFp \ ( T 



(30) 



where the a's and /3's are corrections of order unity 
for which we use (Yakovlev & Levenfish 1995) 



0.59 
: 1.06 



Pn 



np 



0.56 
0.66 



■pp 



0.11 /3pp = 0.70. 



It is important to note that the emissivities of 
the modified Urea and bremsstrahlung processes 
have not been accurately calculated, particularly 
at the highest densities reached in the center of 
the stars we are considering. Voskresensky & Sen- 
atorov (1986) have proposed that when the den- 
sity approaches the critical density for the onset 
of charged pion condensation, the softening of the 
pion mode induces a strong increase in the above 
emissivities. Since this approach assumes the oc- 
currence of charged pion condensation, we will 
not consider it here as part of the minimal sce- 
nario. Nevertheless, less dramatic medium effects 
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arc certainly at work. Recently, Hanhart, Phillips 
& Reddy (2001); Van Dalen, Dieperink & Tjon 
(2003); Schwenk, Jaikumar fc Gale (2004) have 
revisited the bremsstrahlung processes, including 
hadronic polarization up to the two loop level, and 
found a reduction of the rates by a factor of about 
4 at saturation density. In view of this, we will at 
first take equations (23), (25), (28), (29), and (30), 
with the quoted a's and /3's at face value, but will, 
in addition, consider the effects of "cranking up" 
or "down" all modified Urea and bremsstrahlung 
emissivities by a significant factor in § 5.6. 

Once the temperature T reaches the pairing 
critical temperature of either the neutrons or 
protons in a given layer of the star, the correspond- 
ing neutrino emission process becomes suppressed 
by the development of an energy gap A(r) in 
the single particle excitation spectrum (see equa- 
tion (16)). Similarly to what happens for the 
specific heat, the neutrino emissivities are sup- 
pressed by factors which vary approximately like 
exp{—l\{T)/kBT). In our calculations, we employ 
the accurate calculations of these various various 
control functions including pre-exponential factors 
(see Yakovlev & Levenfish (1995) for details) for 
the two modified Urea processes and the three 
bremsstrahlung processes, in the presence of n 
or pairing and/or p ^Sq pairing. Two repre- 
sentative examples, for the neutron branch of the 
modified Urea process, are plotted in the central 
panel of Figure 13: for this specific case neutron 
pairing has a much stronger effect than proton 
pairing since three neutrons, but only one proton, 
participate in the reaction. 

As the temperature begins to approach Tc, new 
channels for neutrino emission through the con- 
tinuous formation and breaking of Cooper pairs 
(Flowers, Ruderman & Sutherland 1976; Voskre- 
sensky & Senatorov 1987) begin to become opera- 
tive. We take emissivities from these "Pair Break- 
ing and Formation" or PBF processes as (Yakovlev 
& Levenfish 1995) 

gP^So = 2.6 X 10^1 f^V^' 

\noJ \mpj 

^^so(w)(^y (31) 



for p^Sq pairing, 

= 1.0 X 10^^ f !^y^' f^'lx 

F^s,iT/n) (y^) (32) 
for n ^Sq pairing, and 

gj:^^^ = 8.6 X 10^1 f 

\nQ J \mnj 

^3p,(T/re)(^y (33) 

for n pairing, assuming again that pairing in 

this last case occurs in the mj = phase. The 
control functions Fis^iT/Tc) and F3p^{T/Tc) are 
shown in the lower panel of Figure 13 and give 
the dependence on Tc- These PBF processes can 
be regarded as nn or pp bremsstrahlung processes 
with a strong correlation in the initial state in the 
case of the breaking of a Cooper pair, or in the final 
state in the case of the formation of a Cooper pair, 
and exemplify an extreme case of medium correc- 
tion to the bremsstrahlung processes. Their effi- 
ciencies are similar to those of the bremsstrahlung 
processes of equations (28), (29), and (30). How- 
ever, they are less sensitive to the values of the 
nucleon effective masses because they are propor- 
tional to m* instead of m**. Furthermore, the T"^ 
dependence of the PBF processes, compared to the 
dependence of the bremsstrahlung processes, 
allows the PBF processes to eventually dominate 
the total neutrino luminosity. The precise value 
of Tc for the neutron and/or proton pairing deter- 
mines when this dominance occurs, as will become 
apparent in § 5.4. The form of the control func- 
tions (Figure 13) show clearly that the PBF pro- 
cess turns on when T reaches Tc, increases its effi- 
ciency as T decreases, and becomes exponentially 
suppressed when the gap approaches its maximum 
size A(0) when T ~ 0.2Tc. 

4. THE NEUTRON STAR ENVELOPE 

It is customary to separate cooling models into 
the interior and the envelope, the latter being the 
upper layer in which a strong temperature gradi- 
ent exists whereas the interior designates every- 
thing inside which becomes isothermal within a 
few years after the birth of the neutron star. 
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Precisely, the envelope can be defined as the 
layer extending from the photosphere, the upper- 
most layer where the emitted spectrum is deter- 
mined, down to a boundary density pi, such that 
the luminosity in the envelope is equal to the total 
surface luminosity of the star, L{r) = L{R). The 
thermal relaxation time scale of the envelope is 
much shorter than the cooling time scale of the in- 
terior so that it can be treated separately as a layer 
constantly in a stationary state. Equation (B2) 
therefore implies that the neutrino emission is neg- 
ligible in the envelope. Since the thickness of the 
envelope is of the order one hundred meters or less, 
the envelope can be treated in the plane paral- 
lel approximation. Within these approximations, 
integration of the heat transport and hydrostatic 
equilibrium equations gives a relationship between 
the temperature at the bottom of the envelope, Tf,, 
and the flux F going through it, or, equivalently, 
a relationship between the effective temperature 
Te and F: F = assTf. This relationship is com- 
monly termed as the "T,, — Ty, relationship" . 

Detailed numerical calculations along this line 
were presented by Gundmundsson, Pethick & Ep- 
stein (1982, 1983) and an analytical approxima- 
tion to these results was provided by Hernquist & 
Applegate (1984), who assumed that the chemi- 
cal composition of matter corresponds to that in 
beta-equilibrium. Gundmundsson, Pethick & Ep- 
stein (1982) found a simple analytical relationship 

Te = 0.87 X 10^ [g, ,,)i/4 (Tb/10«K)0-5'^ , (34) 

where is surface gravity gg measured in 10^^ 
cgs units. (As a rule of thumb, this gives T^, oc 
tI''^ and Te ~ 10<5 K when n ~ 10^ K.) This 
equation illustrates the fact that the dependence 
of the envelope structure on M and R is entirely 
contained in Qg and that Te/g\^^ is independent of 
M and R. This allows us to use "generic" envelope 
models and glue them to the upper layer of any 
stellar model. 

The Sensitivity Strip: Effects of Chemical 
Composition and Magnetic Fields 

The most important finding of Gundmundsson, 
Pethick & Epstein (1982, 1983) is that the T, - 
T(, relationship is mostly determined by the value 
of the thermal conductivity A in a thin layer in 
which ions are in the liquid phase and where A 



is dominated by electron conduction. This layer 
was thus called the "sensitivity strip" in the p — T 
plane. The sensitivity strip is located at lower 
densities for lower temperatures and spans about 
one and a half order of magnitude in density depth. 

The presence of light elements (e.g., H, He, C 
or O) in the envelope can significantly affect the 
Te — Tj relationship if the sensitivity strip is pop- 
ulated by these elements (Chabrier, Potekhin & 
Yakovlev 1997). Lighter elements will burn into 
heavier ones in the thermonuclear regime at high 
enough T and in the pycnonuclear regime at high 
enough p. but conditions in the envelope are usu- 
ally such that H may be present up to densities 
~ 10^ gm cm~^, He up to ^ 10® gm cm~^ and 
C up to ~ 10^^° gm cm~^. The critical tempera- 
tures for thermonuclear burning and densities for 
pycnonuclear burning are well within the sensi- 
tivity strip and one can thus expect a strong ef- 
fect of light element presence on the T;, — Te re- 
lationship. This problem was studied in detail 
by Chabrier, Potekhin & Yakovlev (1997) and 
Potekhin, Chabrier & Yakovlev (1997), who per- 
formed numerical calculations of envelope struc- 
ture with a pure iron plus catalyzed matter chem- 
ical composition and with various amounts of light 
elements. These authors found that the presence 
of light elements can significantly raise the surface 
temperature Tg for a given if they are present in 
sufficient amounts. The larger the amount of light 
elements present, the higher the temperature at 
which their effect will be felt due to the tempera- 
ture dependence of the location of the sensitivity 
strip. But at very high temperatures, the light ele- 
ments have practically no effect because they can- 
not penetrate deep enough. The resulting — 
relationships for various amounts of light elements 
arc shown in Figure 14. 

The presence of a magnetic field can also af- 
fect the structure of the envelope (Greenstein & 
Hartke 1983). The effect is to enhance heat trans- 
port along the field and inhibit transport along 
directions perpendicular to the field. This results 
in a nonuniform surface temperature distribution, 
with a very cold region in which the field is al- 
most tangential to the surface as, e.g., around the 
magnetic equator for a dipolar field, and a corre- 
sponding modification of the — Tt relationship 
(Page 1995). However, the overall effect is not 
very large, but is somewhat sensitive to the pres- 
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ence of strongly nondipolar surface fields (Page & 
Sarmiento 1996). For a field of the order of 10^^ - 
10^^ G, one obtains a slight reduction of com- 
pared to the field-free case, whereas for a higher 
field begins to be enhanced. The enhancement 
of is, however, much smaller than what is ob- 
tained by the presence of light elements (Potekhin 
et al. 2003). Moreover, there are possible instabil- 
ities due to the non-uniformity of the temperature 
(Urpin 2004) which have not yet been taken into 
account in magnetized envelope calculations and 
may somewhat affect these results, but we do not 
expect significant changes. Hence, the important 
case for our purpose would be the maximal re- 
duction of Tg°° obtained for a pure heavy element 
envelope at B = 10^^ G, which is illustrated in 
Figure 14. 

One must finally mention that our calculations 
are based on the assumption of spherical symme- 
try in the interior and that the only asymmetries 
considered, due to the presence of a magnetic field, 
are within the envelope and hence included into 
this outer boundary condition. However, this as- 
sumption is questionable in some magnetic field 
configurations where the field is confined to the 
stellar crust. As shown by Geppert, Kuecker & 
Page (2004), the crust is highly non-isothermal in 
such cases and this can affect the thermal evo- 
lution because the resulting photon luminosity is 
lowered compared to the isothermal crust case. 

5. A GENERAL STUDY OF NEUTRON 
STAR COOLING WITHIN THE "MIN- 
IMAL SCENARIO" 

In this section, we will consider the individual 
effects of the chief physical ingredients which en- 
ter into the modeling of the cooling of an isolated 
neutron star. Our purpose here is twofold: 

(1) to determine the sensitivity of results to uncer- 
tainties in input physics in order to obtain a broad 
range of predictions which, we hope, encompasses 
all possible variations within the minimal cooling 
scenario; and, 

(2) to provide us with the means to identify the 
types of models that will result in the coldest pos- 
sible neutron stars within this paradigm. 

Theoretical refutations of the critical physical 



ingredients needed for these coldest models could 
allow us to raise the temperature predictions and 
possibly provide more, or stronger, evidence for 
"enhanced cooling" . The task of identifying the 
minimally cooling coldest star will be taken up in 
$ 6. An object colder than such a star could be 
considered as evidence for the presence of physics 
beyond the minimal paradigm. 

All results in this section use stars built using 
the APR EOS, except for § 5.8 where the effects of 
the EOS are studied for a star of 1.4 Mq and for 
§ 5.7 where eff'ects of the stellar mass are studied. 

5.1. Neutrino vs. Photon Cooling Eras 
cind the Effect of the Envelope 

The basic features of the thermal evolution of 

a neutron star can be easily understood by con- 
sidering the global thermal energy balance of the 
star 



dE, 



dT 



(35) 



dt dt 

where ii^th is the total thermal energy content of 
the star and Cy its total specific heat. This equa- 
tion is accurate when the star is isothermal, which 
is the case for ages larger than a few decades. Since 
the dominant neutrino processes all have a tem- 
perature dependence, the neutrino luminosity can 
be expressed as 



(36) 



Furthermore, most of the specific heat comes from 
the degenerate fermions in the core for which 



Cv = CT 



(37) 



in the absence of pairing interactions. The photon 
luminosity can be written as 



(38) 



where Te, the effective temperature, is converted 
into the internal temperature T through an enve- 
lope model with a power-law dependence: Te oc 
rpo.b+a ^j^j^ a < 1 (see equation (34) and Fig- 
ure 14 ). Equation (35) is easily integrated in the 
dominantly neutrino and photon cooling eras. 
(1) The neutrino cooling era (L^ ^ L^): In this 
case, 



C 



1 



(39) 
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where To is the initial temperature at time to = 0. 
For T < To, this gives 



6N 



and Tp (xt 12 . 



(40) 



The very small exponent in the evolution dur- 
ing neutrino cooling is a direct consequence of the 
strong temperature dependence of L^,. 
(2) The photon cooling era {L^ ^ L^,): In this 
case, 



* = *i + (7^ - ] (41) 



where Ti is the temperature at time t\. When 
f>ti and T < Ti, we have 



T = 



C 
AaS 



t~^ and Te oc . (42) 



Since a <C 1, we see that, during the photon coo- 
ing era, the evolution is very sensitive to the na- 
ture of the envelope, i.e., a and S, and to changes 
in the specific heat, as induced by nucleon pairing. 

Figure 15 shows the evolution of T^, T^^^^^., 
L^, and L'^ in a simplified model in which no 
pairing has been included, but two extreme cases 
of envelope chemical composition, iron-like ele- 
ments and light elements, arc considered. The 
curves of panel C are analogous to the curves 
of panel A, since L~ = AnR'^asBT^'^. For both 
envelope models the vs. t and 7^^,^^^,^. vs t 
curves follow power laws (i.e., straight lines on a 
log — log plot) in both the neutrino cooling and 
photon cooling eras. For t ~ 10'' yrs, both mod- 
els have the same T^^^^g^ because the envelopes do 
not contribute to energy loss, neither through neu- 
trino emission (due to their low density and very 
small mass) nor through photon emission (since 
<C -^^)- At these times, the model with a 
light element envelope, however, has a higher Tg°°, 
and thus L^, due to the more efficient transport 
of heat in this envelope and will consequently shift 
from neutrino to photon cooling at a much earlier 
time. This trend will not be modified by the in- 
clusion of more realistic physics. During the neu- 
trino cooling era, Tg°° simply follows the evolution 
of the interior temperature and models with light 
element envelopes appear hotter to an external ob- 
server than models with an iron-like envelope, but 



they enter the photon cooling era sooner and sub- 
sequently cool much faster. Neutron stars with 
lesser amounts of light elements in the envelope 
will cool intermediately between the extremes of 
light-element and heavy-element dominated atmo- 
spheres, as displayed in Figure 16. 

5.2. Time Evolving Envelopes 

We consider here the possibility of time evolu- 
tion of the chemical composition of the envelope. 

We assume that the mass of the envelope consist- 
ing of light elements decays with time as 



AML(t) = e-*/^AML(0) 



(43) 



where AMl(O) is the initial mass in light elements. 
This decay could be due to the pulsar mecha- 
nism which injects light elements into the magne- 
tosphere or due to nuclear reactions which convert 
these elements into heavy ones (Chang & Bildsten 
2003a, b). One can expect from this that the star 
will shift from a cooling trajectory corresponding 
to a light element envelope toward a trajectory 
with heavy elements envelope as AA^l decreases. 
Figure 17 illustrates this evolution and shows that 
this shift happens in a short time in the case of an 
exponential mass reduction. This fast evolution is 
in agreement with the Tg — Tb relationships shown 
in Figure 14 where one sees Tg changing from the 
light element case to the heavy element one within 
a small range of variation of AMl, at a value of 
Th depending strongly on AMl. 

5.3. Basic Effects of Pairing 

In this subsection, we briefly illustrate the first 
two significant effects of pairing, suppressions of 
Qi, and Cy. The third effect, neutrino emission 
through the PBF process is studied in the next 
section. The continuous lines in Figure 18 com- 
pare thermal evolutions of the same neutron star 
with and without nucleon pairing (the gaps have 
been chosen so as to maximize effects of suppres- 
sion). The results are very natural: during the 
neutrino cooling era the paired star cools more 
slowly than the unpaired one since its neutrino lu- 
minosity is severely suppressed whereas during 
the photon cooling era it cools faster due to its 
much reduced specific heat. During the neutrino 
cooling era, the suppression of is present, but 
its effect is not so dramatic for three reasons: 
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(1) the lepton contribution to Cv is not sup- 
pressed, whereas is reduced by many orders 
of magnitude as only the very inefficient electron- 
ion bremsstrahlung process in the crust is not 
suppressed, 

(2) when T is not very much less than T^, as is 
partially the case in this example during the neu- 
trino cooling era, the suppression of is preceded 
by a phase of enhancement (see, e.g., Figure 12), 
and 

(3) the cooling curve has a relatively small slope 
when T ~ Tp. From equation (39), and as repre- 
sented schematically in Figure 19, one sees that 
the shift to the T oc t~^/^ power law occurs at a 
time <o-!/ determined by the ratio C/6N. This ra- 
tio is increased by pairing and results in a delayed 
shift, but this amounts to an horizontal translation 
of the cooling curve and hence shows no spectacu- 
lar effect. In contradistinction, during the photon 
cooling era, the shift in the transition time t^-'y 
from neutrino to photon cooling (this occurs ear- 
lier with pairing than without pairing due to the 
smaller value of C/AaS in equation (40)) acts on 
a power law evolution with large slope. 

However, the neutrino emission by the PBF 
process (artificially turned off in this example) al- 
ters significantly these simple results as can be 
seen from the dotted line in Figure 18, and is de- 
scribed below. 

5.4. The PBF Neutrino Process 

We consider here in detail the effect of the PBF 

neutrino process. Given the strong T^F(T/A) 
temperature dependence and the density depen- 
dence of A, the overall effect can only be assessed 
by complete calculations presented here and in § 7. 

As a first step, we consider separately the tem- 
perature dependence of the luminosities due to the 
n and p^S^ gaps in the core of a 1.4 Mq star 
built with the EOS of APR. Results for four differ- 
ent n gaps are shown in Figure 20. The lower 
panel shows the Tc profiles of these four gaps as 
a function of the volume of the core (left hand 
scale) and the density (right hand scale). A ver- 
tical line in this panel, which corresponds to an 
isothermal core, allows one to visualize the amount 



of the core's volume which is paired. The upper 
panel shows the corresponding PBF neutrino lu- 
minosity L^^^. Also plotted are the surface pho- 
ton luminosity corresponding to an iron-like en- 
velope (dotted line) and the total neutrino lumi- 
nosity from the modified Urea and bremsstrahlung 
processes without pairing suppression. Notice that 
when T ~ 10^ K, the star is usually not isother- 
mal: the crust is warmer than the core and thus 
is larger than indicated in this figure. 
When T decreases, L^^^ grows very sharply 
for each gap when T reaches the maximum Tc of 
the gap reached in the core (PBF neutrino emis- 
sion turns on) and then decreases with a T depen- 
dence which is between a and a power law. 
This results from the overall T'^F{T/A) temper- 
ature dependence of the PBF neutrino emissivity 
combined with the density dependence of Tc which 
determines how much of the core, at this given T, 
is contributing to L^^^ . Once T is much below 
the minimum value of Tc reached in the core, L^^^ 
becomes exponentially suppressed. In the cases 
of gaps "b" and "c" , this suppression appears at 
T < l{f K, for gap "a" when T <C 2 x 10^ K, 
whereas for the gaps "T72" this suppression does 
not appear since Tc reaches very low values and 
for any T there is always a significant volume of 
the core where T r^Tc- 

When T - 10** K, L.^ dominates over so that 
the important range of T to consider is 10® — 10^ 
K and in this range the figure shows clearly that 
the relatively small n gaps as "T73" and 
"a" generate a L^^^ which is higher, by about 
one order of magnitude, than the combined 
that the modified Urea and bremsstrahlung pro- 
cesses would produce when no pairing is present 
(dashed line in the figure). The modified Urea and 
bremsstrahlung processes are of course strongly 
suppressed in the presence of pairing, but these re- 
sults show that, compared to the no pairing case, 
pairing can actually increase the total neutrino lu- 
minosity through the PBF neutrino emission if the 
gaps have the appropriate size. 

Very similar results are obtained when consid- 
ering the PBF neutrino emission due to the p^Sq 
gaps as shown in Figure 21. The three differences 
with respect to the n case are, first, that no 
calculation of the p ^Sq gap reaches a Tc as high 
as our case "c" for the n ^P2, second, all p ^Sq 
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gaps vanish in the inner part of the core which im- 
plies that the suppression of L^^^ at low T does 
not show up and, third, the volume of the core 
in the superconducting state is generally smaller 
than the volume in the n superfluid state, result- 
ing in lower values of L^^^ for p than for n. 

5.5. Comparison of Various Neutrino Lu- 
minosities 

Having compared, in the previous subsection, 

the neutrino luminosities from PBF assuming an 
isothermal interior, we now proceed to analyze 
them, together with other processes, in realistic 
cooling calculations which take into account the 
exact temperature profile within the star. The re- 
sults are shown in Figure 22. We use a 1.4 Mq star 
built with the EOS of APR and we fix the n and 
P^Sq gaps considering that the major uncertainty 
in the neutrino luminosity is due to uncertainty in 
the size of the n gap (see the results of §5.4). 
We consider the three cases "a" , "b" , and "c" for 
the magnitude of the n ^P^ gap. 

The three panels of Figure 22 show c;learly that 
at very early times the cooling is driven by the 
modified Urea and nucleon bremsstrahlung pro- 
cesses, but that once pairing occurs the neutrino 
emission from the PBF process takes over because 
of its efficiency and because the other processes 
are suppressed. 

At ages relevant for the presently available 
data, t ~ 10^ yrs, we find in all cases that the PBF 
neutrinos are the main cooling agent, until photon 
emission takes over sA, t ^ 10^ years. There is a 
competition between the neutrino emission from 
n^P^ and p^Sq pairing: the smallest gap drives 
the cooling between 10^ and 10^ years. With the 
assumed p^So gap from AO, we see that in case 
"a" the n^P2 gap drives the cooling (Figure 22 
left panel), whereas in case "b" or "c" (Figure 22 
central and right panels, respectively) the proton 
pairing drives the cooling. 

This is in agreement with the results of Fig- 
ure 20 which show that the case "a" n ^P2 gap 
is the most efficient and that even smaller gaps, 
as in T72, do not result in a significant enhance- 
ment of the PBF neutrino luminosity or do it too 
late, i.e., at a time when photon emission domi- 
nates the cooling. Comparison of Figure 20 with 
Figure 21 shows that in the case of a larger n^P2 



gap, one can expect that the p ^Sq, PBF neutrinos 
will dominate the cooling, in agreement with our 
findings of this subsection. Finally we compare 
cooling trajectories with our three n '^P2 gaps and 
a vanishing gap in Figure 23 explicitly confirm- 
ing that n ^P2 gaps with Tc of the order of 10^ K 
are the most efficient gaps with regard to neutrino 
cooling through the PBF process. 



23 



Density [10^"^ g cm""^] 
987 65 432 



-3n 



I 



E 
o 



(U 

o 

CM 

o 



o 



Density [10^"^ g cnn 



-3i 




^ 0.5 -P 



Stellar Volume [10^** cm"^] 



Fig. 11. — Cumulative specific heats of e, fj,, p, and 
n vs. stellar volume in the core of a 1.4 Mq star 
built using the APR EOS, at temperature T = 10^ 
K. Nucleons are assumed to be unpaired. 
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Fig. 12. — Specific heat of neutrons in the core 
of a 1.4 Mq star built using the APR EOS, at 
six different temperatures, from 0.99 down to 0.1 
times 10^ K illustrating the effects of pairing. No 
gap is present for results shown by the dotted 
lines. The n gap "a" is assumed for continu- 
ous curves. This gap has a maximum Tg of 1 x 10^ 
K at /9 = 3.61 X 10^^ gm cm~^. the continuous 
curves show the jump of Cy by a factor 2.188 (see 
Figure 13) at the two zones where T = Tc and its 
progressive suppression when T <C Tc in the layers 
in between. 
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Fig. 13. — Control functions for pairing effects in 
^5*0 (for neutrons and/or protons) and (neu- 
trons) channels. The top panel shows the function 
rclavant for the specific heat (§ 3.4), the central 
panel that for the neutron branch of the modified 
Urea process with either p ^Sq ov n^P2 pairing 
(§ 3.5) and the bottom panel is for the PBF pro- 
cess (§ 3.5). 




Fig. 14. — Relationship between the effective 

temperature and the interior temperature T;, 
at the bottom of the envelope assuming various 
amounts of light elements parametrized by 77 = 
g1i^l\Mi^/M (AMl is the mass in light elements 
in the envelope, QsU the surface gravity in units 
of 10^* cm s~^, and M is the star's mass), in the 
absence of a magnetic field (Potckhin, Chabrier 
& Yakovlev 1997). Also shown are the - U 
relationships for an envelope of heavy elements 
with and without the presence of a dipolar field of 
strength of 10^^ G following Potekhin & Yakovlev 
(2001). 
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Fig. 15. — Neutrino and photon cooling eras for 

two models of non magnetized envelopes formed 
by heavy iron-like elements (labeled "H") and a 
maximum amount of light elements (labeled "L" ) . 
The effective temperature (panel A), the cen- 
tral temperature (panel B) and neutrino (contin- 
uous lines) and photon (dotted lines) luminosities 
(panel C) , all redshifted to infinity, are shown as a 
function of time. Pairing effects are not included 
in these calculations. 
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Fig. 16. — Effect on the cooling of various amounts 
AMl (in solar masses) of light elements in the 
envelope. The two dashed curves, "H" and "L", 
are the same as in Figure 15. Pairing effects are 
not included in these calculations. 




Log t [yrs] 

Fig. 17. — Transition of cooling trajectories be- 
tween a model with a heavy element envelope 
(doted curve labeled "H" ) and a light element en- 
velope of maximum mass (dotted curve labeled 
"L"). Continuous curves show evolution of mod- 
els with "decaying" envelopes (see equation (43)) 
with various decay times r as indicated. Also in- 
dicated are masses of the light element envelopes 
at the moment the star begins its shift toward the 
heavy element envelope trajectory. 
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Fig. 18. — Comparison of the cooHng of a 1.4 

M© star, built using the EOS of APR, without 
and with nucleon pairing. In the model with pair- 
ing neutrino emission by the PBF process has been 
either artificially turned off (continuous line) or al- 
lowed (dotted line). Neutron ^5*0 pairing is from 
AWP, pairing from our case "c" and p ^Sq 
pairing from AO, as labeled in Figures 8, 10, and 
9, respectively. The envelope is assumed to be 
composed of heavy elements. 
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Fig. 20. Upper panel: neutrino luminosities vs 
temperature from the PBF process for four dif- 
ferent n gaps labeled as in Figure 10. Lower 
panel: Tc for the four neutron gaps vs density 
and enclosed volume. 
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Fig. 19. — Schematic representation of the power- 
law cooling behaviors and the effect of pairing. 
The time to-i, denotes when the central tempera- 
ture falls to a value small enough that T oc t~^^^ 
becomes valid. The time r,^_-y denotes the transi- 
tion from neutrino to photon cooling eras. 
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Fig. 21. — Upper panel: neutrino luminosities vs 
temperature from the PBF process for 3 different 
p ^Sq gaps labeled as in Figure 8. Lower panel: 
critical temperature Tc for the 3 p gaps vs 
density and enclosed volume. 
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Fig. 22. — Comparison of luminosities from various processes during three realistic cooling histories: photon 
("7")' 3.11 J/-processes ("Total f"), modified Urea and nucleon bremsstrahlung ("MU-Br. ;/"), and PBF 
("PBF u") from n and p ^Sq pairing marked by "n" and "p", respectively. PBF neutrinos from the n 
^Sq gap are not shown explicitly, since their contribution is always dominated by other processes, but they 
are included in the total v luminosity. In all three cases, the p ^Sq gap is from AO, the n ^Sq gap from 
AWP, whereas the n^P2 gap is our model "a" (left panel), "b" (central panel) and "c" (right panel). 
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5.6. More-Modified Urea Cooling 

The modified Urea and bremsstrahlung pro- 
cesses all involve four nucleons and are processes 
in which energy-momentum transfer occurs via 
strong interactions in the medium. The associ- 
ated emissivities are sensitive to one's assumptions 
about in-medium strong interactions and their efS- 
cicncies arc difficult to assess with certainty. Given 
this, we consider it important to study the effect of 
this uncertainty in a simple, but drastic way: we 
simply multiply the gMUrca .^^^^ q^'^°'^ emissivities 
by a constant factor F, taking F to be 1/10, 10 
or 100. A factor 1/10 or 10 could be acceptable, 
whereas a factor 100 is probably exaggerated. 

The results of Figures 20 and 21 showed that in 
the presence of n and p^So pairing, most rea- 
sonable gaps produce a neutrino emission by the 
PBF process which is much more intense than the 
modified Urea and bremsstrahlung in the absence 
of pairing by at least one order of magnitude in 
the important temperature range T 10® — 10^ 
K. Thus, in a realistic calculation a factor F = 10 
is not expected to lead to a significant change 
in the cooling. This is confirmed by our results 
shown in Figure 24. The models with no pairing 
clearly show enhanced cooling when = 10 or 
100, and reduce cooling when F = 1/10, whereas 
when pairing is included the models with F = 10 
and 1/10 are practically indistinguishable from the 
unenhanced case and only the, probably unrealis- 
tic, case F = 100 leads to a faster cooling. 

These results are important and fortunate, 
since they show that the uncertainty in the ac- 
tual efhciency of the modified Urea rate has no 
significant effect on the predictions of the minimal 
scenario when pairing, and the corresponding neu- 
trino emission from the PBF process, is included 
in a realistic way. 

5.7. Effects of Neutron Star Mass 

In the case that neutrino cooling occurs only 
through the modified Urea and bremsstrahlung 
processes, as required by the tenets of the minimal 
cooling scenario, one can expect that the cooling 
curves in the neutrino cooling era will show prac- 
tically no variation with neutron star mass, be- 
cause there are no energy or density thresholds 
for these processes. This situation will change 
drastically for the case in which enhanced cool- 
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Fig. 23. — Comparison of cooling trajectories with 
vanishing n gaps, labeled "0", and our three 
model gaps "a", "b", and "e" (see Figure 10). The 
n^So gap is from AWPIII and the p^So gap from 
AO (see Figures 8 and 9). Results are for a 1.4 
Mq star built using the EOS of APR with a heavy 
element envelope. 




Fig. 24. — Cooling with adjusted modified Urea, 
for F = 1/10, 1, 10 and 100 as indicated, with and 
without nuclcon pairing. Assumed pairing gaps 
are from our model "a" for n ^P^ and from AO 
for p^Sq (and n^So pairing from AWP for which 
the effect is very small). The envelope is assumed 
to be composed of heavy elements. 
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ing through direct Urea processes becomes possi- 
ble either through nucleons or due to the presence 
of exotica. 

Figure 25 confirms that in the absence of pair- 
ing, there is almost no mass effect, both during 
the neutrino and the photon cooling era. Simi- 
larly, when n, but not p, pairing is included, the 
mass dependence is also small, though larger than 
with no pairing at all. When the p gap is included, 
the main variation with mass occurs in the pho- 
ton cooling era in which more massive stars cool 
more slowly. This is a direct consequence of the 
lesser suppression of the proton specific heat with 
increasing mass, since the p ^Sq gap vanishes at 
high density and there is an increasingly larger un- 
paired region when M increases (Cvip) is larger 
for larger M). The chosen n gap reaches the 
center of the star in all cases and thus Cv{n) is 
strongly reduced for all masses, which explains the 
small mass dependence when only n gaps are taken 
into account. In the case that the n^P2 gap would 
also vanish at high density, we would obtain an ad- 
ditional mass dependence. 

5.8. Effects of the Equation of State 

In exploring the high density EOS, one can ex- 
pect three sources of effects: 

(1) general relativistic effects due to change in the 
star's compactness, 

(2) differences in the n and p effective masses, and 

(3) differences in the volume of the star in the 
various paired states. 

Figure 26 shows results for the four EOS's se- 
lected in § 3.1. When no pairing is included, there 
is little variation with the EOS, and slight varia- 
tions exist when pairing is considered. The rea- 
sons are essentially the same as those discussed in 
conjunction with the stellar mass (see the previous 
subsection) and are due to the density dependence 
of the p^So gap and, to a much lesser degree, that 
of the n gap. The very small differences in 
the unpaired models simply reflect that the four 
chosen EOS's are rather similar because of con- 
straints imposed by the minimal cooling scenario: 
the differences in the stars' compactness and nu- 
cleon effective masses are very small. 




Log t [yrs] 



Fig. 25. — Effects of the stellar mass: cooling 
of stars of various masses built using the EOS of 
APR, with and without pairing. Models with pair- 
ing have u^Sq gap from AWP and n^P2 gaps from 
our model "a' and either no n '^P2 gap (dashed 
curves) or p gap from AO (continuous curves). 
Stellar masses are indicated in the cases with the 
three types of pairing, whereas for similar cases 
without proton pairing or with no pairing at all 
the trajectories are too similar to be separately 
labeled. The envelope is assumed to be composed 
of heavy elements. 
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6. MINIMALLY COOLING COLDEST 
NEUTRON STARS 

One of the main goals in this work is to deter- 
mine how cold an observed neutron star should 
become to be incompatible with the predictions 
of the minimal scenario. Armed with the results 
of the previous section, we can now identify the 
fastest cooling models within this scenario. 

6.1. Neutrino cooling era 

During the neutrino cooling era, Figure 23 
shows that the lowest Te's are obtained due to the 
PBF process when the n gaps are of the size 
of model "a", i.e., with T^s of order at most 10^ 
K in most of the stellar core. The p^Sq gaps can- 
not compete with the most efficient n gaps, 
because proton gaps are restricted to a smaller 
volume; compare Figures 20, 21, and 22. These 
fastest neutrino cooling models have a very weak 
dependence on the mass of the star (Figure 25). 
These models also require that the envelope be 
made of heavy elements or, if it contains light ele- 
ments, their amounts should be much smaller than 
10-11 Mq (see Figure 16). 

6.2. Photon cooling era 

The physical processes that control cooling in 
the photon cooling era arc quite different from 
those in the neutrino cooling era. Neutrino emis- 
sion from any of the possible processes make only 
a small contribution in the photon cooling era. 

The two crucial ingredients are the envelope, 
which determines the photon luminosity and the 
specific heat, which controls pairing (see $ 5.1). 
A light element envelope, producing a higher Tg 
and hence a higher L for a given core tempera- 
ture, leads to fast cooling; an amount above 10"^ 
Mq of these elements is necessary (see Figure 16). 
Concerning the total specific heat, the strongest 
reduction can be achieved by pushing baryon pair- 
ing to the extreme, and this means considering low 
mass neutron stars so that the p^Sq gap is more 
likely to reach the center of the star. Pursuing the 
trend indicated in Figure 25, we consider effects 
of the various p ^Sq gaps of Figure 9 for a low 
mass, 1.1 Mq , neutron star. Results are shown 
in Figure 27. The proton kp at the center of this 
star has a value of 1.1 {m~^ (see Figures 5 and 6); 
the inset of Figure 27 shows a direct mapping of 



the density at which the p ^Sq gap vanishes with 
Te at these times. The fastest cooling model cor- 
responds to the p^Sq gap "CCDK," which has a 
Tc of 1.44 X 10^ K at the center of the star and 
hence produces a complete suppression of the pro- 
ton specific heat in the photon cooling era. A-P^Sq 
gap with a higher Tc at the center of the star, or a 
gap that vanishes at higher densities (not reached 
in this star), would lead to the same cooling curve. 
Similar considerations apply to the n ■^P2 gap. To 
illustrate this, we used our model gap "a" in Fig- 
ure 27. Any other gap with a Tc higher than a few 
times 10* K would result in the same total suppres- 
sion of the neutron specific heat and, therefore, to 
exactly the same cooling curve. 

These results will be important when compar- 
ing our predictions with data in the next section, 
particularly for young stars with ages of the order 
of a few times 10'' years, such as the Vela pulsar 
and PSR 1706-44. 

7. COMPARISON OF THE MINIMAL 
COOLING SCENARIO WITH DATA 

In § 5, we analyzed in some detail the effect of 
each physical ingredient that shapes the cooling 
history of a neutron star within the minimal sce- 
nario. In § 6, we identified the fastest cooling neu- 
tron star models in this scenario. The combined 
effects of these ingredients in realistic models, to- 
gether with comparisons to the presently avail- 
able temperature and luminosity measurements, 
are presented below. 

Our task is greatly simplified by the fact that 
the EOS is considerably constrained by the tenets 
of the minimal scenario (see results of § 5.8). 
Moreover, as shown in § 5.7, the precise mass of 
the neutron star also has little effect, with the pos- 
sible exception of low mass stars at ages around a 
few times 10** years. (This, of course, is changed 
drastically once we go beyond the minimal sce- 
nario and allow for enhanced neutrino emission 
processes to occur at high density.) We can there- 
fore restrict our attention mostly to the thermal 
evolution of a "canonical" 1.4 Mq neutron star 
built with the EOS of APR. In contrast, the chem- 
ical composition of the envelope and the extent of 
pairing of both neutrons and protons will play sig- 
nificant roles. 

As shown in § 5.1, the presence of light elements 
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Fig. 26.— Effects of the EOS: cooling of 1.4 
Mq stars built using the four chosen EOS's. 
EOS's are labeled in the cases with pairing, 
whereas for similar cases without pairing the tra- 
jectories are too similar to be separately labeled. 
Pairing gaps: n ^Sq from AWP, n from our 
model "a" and p ^Sq from AO. The envelope is 
assumed to be composed of heavy elements. 
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Fig. 27. — Cooling of a 1.1 Mq star (continuous 
lines) and with various p ^So gaps (as labeled in 
the inset): 1 - NS, 2 - T, 3 - AO, 4 - BCLL, 5 
- CCY_ms, and 6 - CCDK (see Figure 9 for no- 
tation). The n ^So gap is from AWP, and the 
n gap from our model "a" . The dotted line 
is the same 1.1 Mq star, but without any pairing 
and the dashed line is for a 1.4 Mq star, with the 
same pairing gaps as in Figure 25. Envelopes are 
all assumed to be composed of light elements. 



in the envelopes of young neutron stars leads to 
effective temperatures that are larger than those 
without any light elements during the neutrino 
cooling era, whereas it implies a faster cooling dur- 
ing the later photon cooling era. Thus, for an 
assumed high density structure of the star, there 
exists a whole family of models limited by the two 
extreme cases of envelopes: those with only heavy 
elements, and those with a maximum amount of 
light elements. Stars with envelopes containing 
only a small amount of light elements will evolve 
on intermediate tracks shifting from a track close 
to the former one toward the latter one as illus- 
trated in Figure 16. Conversely, stars can evolve 
in the opposite direction if the envelope composi- 
tion changes with time from light elements toward 
heavy elements; such an evolution is very abrupt 
as shown in Figure 17. 

On the other hand, the occurrence of pairing 
also accelerates cooling during the photon cooling 
era through the reduction of the specific heat (see 
§ 5.3), whereas pairing effects during the neutrino 
cooling era are more delicate. Neutrino emission 
from the modified Urea and bremsstrahlung pro- 
cesses is suppressed, but the breaking and forma- 
tion of Cooper pairs can easily, with appropriate 
gaps, become vastly more efficient than the former 
processes. As a result, depending on its size and 
density dependence, pairing can lead to faster or 
slower cooling during the neutrino cooling era; its 
effect has to be considered carefully. 

Our main results arc presented in Figure 28 
and compared with data. For the reasons dis- 
cussed in § 2.3, we present them in two forms: 
effective tcmpcrat\.ire T,°° versus time and lumi- 
nosity versus time. We divide our results into 
three subclasses depending on the size of the n 
gap, giv(ui that this gap is the most uncertain 
one: a vanisliingly small gap and our gap mod- 
els "a" and "b" . Figures 22 showed that case "c" 
results in neutrino emissions very similar to that 
of model "b" and Figure 23 confirmed that the 
resulting cooling trajectories are practically iden- 
tical for these two large n gaps and we there- 
fore do not need to include results for the gap "c" 
here. For each assumed neutron gap, it is still 
necessary to consider uncertainties in the n and p 
^So gaps. Varying these gaps is less dramatic than 
varying the n ^P2 gap and we consider 15 differ- 
ent combinations (see the caption of Figure 28) 
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which we plot together. We obtain, for each as- 
sumed n gap, two sets of very closely packed 
curves, one for each envelope composition. The 
size and extent of the n ^Sq gap has very little ef- 
fect, since it is mostly restricted to the crust and 
encompasses only a small part of the star's volume. 
This leads to small differences in the early cooling 
when the star has not yet reached isothermality, 
at ages ~ 3 — 100 years, and the surface tempera- 
ture is still controlled by the evolution of the crust. 
Among the p ^Sq gaps, the ones which can reach 
higher densities will lead to slightly faster cool- 
ing both during the neutrino cooling era, because 
of the enhanced neutrino emission from the PBF 
process, and during the photon cooling era, due to 
the resulting smaller specific heat. 

Figure 23 demonstrated that models with the 
n^P2 gap "a" yield the coldest stars and Figure 28 
shows that the spread in results due to the vari- 
ation of the other two gaps is much smaller 
that in the other two scenarios. This is because the 
n^P2 gap "a" maximizes neutrino emission by the 
PBF process (see Figure 20) and because the neu- 
trino luminosity due to the proton PBF process 
is lower (compare Figure 21 with Figure 20). We 
have studied many other models with a n gap 
similar to our case "a" , but with slightly different 
maximum values of Tc and different density depen- 
dences, and discerned negligible differences. We 
are thus confident that the results presented here 
reflect the smallest temperatures possible within 
the constraints of the minimal scenario. We can 
obtain a slightly faster cooling in the photon cool- 
ing era for low mass neutron stars (near 1 to 1.2 
Mq) as discussed in § 6, and this case will be pre- 
sented separately at the end of this section. 

We now compare observational data for specific 
neutron stars with the suite of models encompass- 
ing the minimal cooling scenario. 

7.1. RX J0822-4247 and IE 1207-52 

These two are young and are the hottest known 
stars. Their inferred temperatures are higher 
than the predictions of all our models with heavy 
element envelopes, but are compatible with all 
models with light element envelopes. This may 
be considered as possible evidence for the pres- 
ence of a significant amount of light elements in 
the upper layers of these stars. However, when 
considering luminosities, for select values of the 



n gap, IE 1207-52 is compatible with hav- 
ing an heavy elements envelope. RX J0822-4247, 
however, remains more luminous than any of the 
heavy-clement envelope models but is compatible 
with light-element envelope models. 

7.2. PSR 0656+14, PSR 1055-52, Geminga, 
RX J1856.5-3754, and RX J0720.4- 
3125 

These are the five oldest observed stars. Fits of 
their spectrum to light-element atmospheres result 
in radii too large to be compatible with the neu- 
tron star hypothesis. Blackbody spectral fits re- 
sult in too-small radii, but it is possible for heavy- 
element atmospheres or two-temperature black 
bodies to be constructed that produce compat- 
ible radii. For consistency, we have restricted 
the data appearing in these plots to be a result 
of either light-element atmosphere fits or single- 
temperature blackbody fits. Inferred tempera- 
tures are more sensitive to the assumed atmo- 
spheric composition than are luminosities, and for 
these five objects, the L versus age plots are prob- 
ably more reliable and representative of the obser- 
vational uncertainties. Consequently, blackbody 
fits result in relative positions for temperatures 
that are quite different than those of the lumi- 
nosities. 

Except for RX J1856. 5-3754, the large uncer- 
tainties on the age and the luminosity of these ob- 
jects preclude definite conclusions. If we consider 
the upper limits to their agea and/or luminosities, 
we find them too bright and must invoke the pres- 
ence of some strong heating process. On the other 
hand, considering the lower limits to ages and/or 
luminosities they appear compatible with the min- 
imal scenario independently of assumptions about 
pairing. 

7.3. Vela, and PSR 1706-44 

Very intriguing objects are the pulsars PSR 

0833-45 ("Vela") and PSR 1706-44. Vela has been 
repeatedly proposed as a candidate for enhanced 
cooling or exotic matter, but our results are in- 
conclusive with respect to these claims. For this 
star, the effective area is compatible with emis- 
sion from almost the entire surface of a neutron 
star and both types of plots, T or _L vs age, are 
equivalent. With an n '^^2 gap chosen to maximize 
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neutrino emission from the PBF process, as in our 
case "a" , the discrepancy of Vela with the theo- 
retical prediction is not significant whereas for a 
vanishing n gap it is very large. However, for 
any non- vanishing 71*^^2 gap and a low assumed 
stellar mass (see Figure 29), several of the light 
element envelope models reach the temperature 
of Vela at an age of 20,000 yrs, i.e., well within 
its age uncertainty and even less than some esti- 
mates of the supernova remnant age, 1.8±0.9x 10^ 
yrs (Aschenbach, Egger & Triimper 1995). These 
models correspond to p ^Sq gaps which extend to 
relatively high densities and hence result in strong 
suppression of the proton specific heat in most of, 
if not all, the core (see $ 6). An interesting fea- 
ture of these models is the very fast decrease of 
temperature 



AT, 



At 

-0.85 — . 

t 



(44) 



which for an age of 20,000 yrs gives a decrease of 
the observable x-ray flux of 0.17% every decade. 
In this case one could interpret these results as 
indicating that the Vela pulsar is a low mass neu- 
tron star with a thick light element envelope and 
in which neutrons and protons are paired in the 
entire core. On the other hand, they could favor 
a neutron star whose core neutrons have a Tc of 
the order of 10*^ K, without any constraint about 
the proton pairing and the stellar mass, but this 
star must have a heavy clement envelope with at 
most 10^^'^ Mq of light elements at the surface 
(see Figure 16). 

PSR 1706-44 is in a similar situation, but with 
larger uncertainties both in T and in L^o, due to 
the large distance and age uncertainties. Confir- 
mation of its association with the supernova rem- 
nant G 343.1-2.3 (McAdam, Osborne & Parkinson 
1993) would help constrain both its distance and 
age. 

7.4. Barely detected and undetected ob- 
jects 

Sources that have negligible or no observed 
thermal emissions, listed in Table 3 and plotted 
in Figure 2, are compared with our results in Fig- 
ure 30. The upper luminosity limits for the two 
objects CXO J232327.8-f-584842 (in Gas A) and 
PSR J0154-I-61 are well within the prediction of 
the minimal scenario whereas the limit of PSR 



J1124-5916 is on the lower side, but still compati- 
ble. 

Most interesting are the two stars PSR J0205-h6449 
(in 3C58) and RX J0007.0-h7302 (in CTA 1) whose 
upper limits are clearly below our predictions. The 
remaining four objects, with no point-like emis- 
sions of any kind observed to date, would provide 
definitive evidence for enhanced cooling if it could 
be shown that neutron stars in fact exist in any of 
them. 
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Fig. 28. — Comparison of predictions of the minimal cooling scenario with data, all models being 1.4 Mq stars 
built using the EOS of APR. Left panels: effective temperature vs. age. Right panels: luminosity Loo 
vs age. The upper, middle, and lower panels correspond to throe different assumptions about the size of the 
n gap as indicated in the panels. In each panel, the two sets of curves correspond to the two extreme 
models of envelope chemical composition: light elements or heavy elements, as labeled in the upper left 
panel. For each set of curves, the 15 different curves correspond to three different choices of the n ^Sq gap 
("AWPII", "AWPIII", and "SCLBL" as labeled in Figure 8) and five different p ^Sq gaps ("CCYms", "T", 
"NS", "AO", and "BOLL" as labeled in Figure 9). 
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8. COMPARISON WITH OTHER WORKS 

The literature on neutron star cooling is ex- 
tensive and dates back to the early 60's. De- 
tailed studies of the "standard" cooling of iso- 
lated neutron stars have been presented, e.g., by 
Nomoto & Tsuruta (1987), van Riper (1988), and 
Schaab et al. (1996), but none of these works 
had included the neutrino emission from the PBF 
process. Models incorporating the PBF pro- 
cess were first presented by Schaab et al. (1997) 
and subsequently by Page (1998). These two 
works, however, did not explore effects of various 
EOS's and/or various sets of pairing gaps within 
the framework of the "standard" scenario. Our 
present results are in agreement with the results 
of Nomoto & Tsuruta (1987), van Riper (1988) 
and Page (1998) when similar input physics are 
used and allow a close comparison. 

We find, however, large differences with the re- 
sults of Schaab et al. (1996) and Schaab et al. 
(1997) during the photon cooling era in which the 
models of these authors cool much more slowly 
than our models. This discrepancy is probably 
attributable to differences in the specific heat. 
For example, during the photon cooling era, the 
Schaab et al. models in which pairing is considered 
cool in almost the same way as models in which 
the effects of pairing are neglected. The models 
without pairing have luminosities about one order 
of magnitude higher than our corresponding mod- 
els at ages in the range 10^ — 10'' years. 

Extensive calculations of the effect of the 
PBF process have been presented by Lcvenfish, 
Shibanov & Yakovlev (1999) in a simplified model 
in which both the p ^Sq and n gaps were 
assumed to be constant in the entire core of the 
star. This simplification of uniform Tg (which im- 
plies that protons are paired in the entire core 
of their stars, whereas in our models the protons 
in the central part of the star are unpaired for 
1.4 Mq stars) has the effect of overestimating the 
neutrino emission by the PBF process. Their re- 
sults are consistent with our findings in that they 
obtain somewhat cooler stars. For example, they 
can reach logT^"" = 5.8 at age t = 10"* years, 
whereas our coolest model with n^P2 gap of our 
model "a" has logT^ = 5.9 at the same age. 

Subsequent works of Yakovlev ct al. (2002) and 
Yakovlev & Haensel (2003), which explored more 



realistic, density-dependent gaps, obtained results 
which are in good agreement with ours. For mod- 
els in which no enhanced neutrino emission is at 
work, these two works obtain a minimal log of 
5.9, as we do, at an age oit = 10^ years. 

Several recent works, Kaminker, Haensel & 
Yakovlev (2001), Yakovlev, Kaminker & Gnedin 
(2001), Kaminker, Yakovlev, & Gnedin (2002), 
and Tsuruta et al. (2002), studied the cooling of 
neutron stars within the "standard" scenario (in- 
cluding the PBF processes), but with enhanced 
neutrino emission at high density. In addition, 
they did not explore the range of parameters that 
are considered here. Consequently, these studies 
cannot be compared to the minimal scenario pre- 
sented here. 

In a recent review, Yakovlev & Pethick (2004) 
have examined standard and enhanced cooling in 
an attempt to fit all the data within one model. 
In this work, the data on individual objects arc fit 
by treating the mass of the star as a free param- 
eter. Models in which the n ^P^ gap is taken to 
be vanishingly small (Schwenk & Friman 2004) are 
favored, together with up ^Sq gap that persists up 
to very high densities and with a Tc of the order 
of 7 X 10*^ K. This latter feature is incompatible 
with all calculations of this gap we report in Fig- 
ure 9, but is not absolutely excluded. Our models 
for 1.1 Mq with nucleon pairing in the entire core 
(see Figure 29) are similar to those of Yakovlev & 
Pethick with non-vanishing n ^P2 gaps, whereas 
OMv models of heavier stars with extensive neutron 
pairing but a large volume of unpaired protons are 
similar to their models with vanishing n ^P2 gaps 
and extensive proton pairing (see Figures 28). 

9. DISCUSSION AND CONCLUSIONS 

We have presented a detailed study of the 

thermal evolution of an isolated neutron star us- 
ing what we term as the minimal cooling sce- 
nario. This scenario is an extension of the well- 
known "standard cooling" scenario to include the 
effects of nucleon pairing and complements neu- 
trino emission by the modified Urea and nucleon 
brcmsstrahlung processes with the pair breaking 
and formation (PBF) process. We have confirmed 
the results of previous works by others that for 
many models of nucleon pairing, the PBF pro- 
cess actually dominates the cooling of the star (see 
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§ 5.4) and hence is an essential ingredient of the 
minimal cooling scenario. 

Among the four parameters we proposed for an 
overall classification of neutron star cooling mod- 
els, we found that the EOS at supranuclear densi- 
ties is well constrained by the requirements of the 
minimal cooling scenario. Moreover, we showed 
that the stellar mass has little effect on the re- 
sults. We emphasize that for scenarios beyond the 
minimal one, i.e., when new particles and neutrino 
emission processes appear, these two parameters 
definitely gain importance. The other two param- 
eters we considered, pairing properties of the nu- 
cleons and chemical composition of the envelope, 
introduce the largest uncertainties in our theoret- 
ical predictions. 

We singled out three subclasses of scenarios due 
to uncertainties in the size and extent of the n 
^P2 gap. For this gap, we considered three dif- 
ferent cases: a vanishingly small gap, a somewhat 
small gap (our model "a") and a relatively large 
gap (model "b"). Within each of these cases, vari- 
ations of the n and p ^Sq gaps covering the pub- 
lished ranges of these gaps were also considered. 

With respect to the chemical composition of 
the envelope, we singled out two extreme cases: 
an envelope consisting of heavy iron-like elements 
and an envelope containing essentially only light 
elements. For each choice of the n gap, we 
obtained two families of closely packed cooling 
curves representing each of the two extreme en- 
velope cases, with some spreads due to variations 
of the n and p ^Sq gaps. Intermediate envelope 
chemical compositions, or its possible temporal 
evolution, result in cooling trajectories interme- 
diate between the extremities (see Figures 16 and 
17). 

Comparing our results with observationally in- 
ferred temperatures Too and luminosities Loo of 
eleven isolated neutron stars (Figures 28 and 29), 
we found that the observations were in overall 
good agreement with the minimal cooling scenario 
taking into account the uncertainties of envelope 
and pairing properties as well as those of the ages 
and inferred Too and Loo for these stars. It is prob- 
ably not possible to understand these data within 
a single model with a unique envelope chemical 
composition. Considering that the compositions 
of the upper layers are strongly dependent on 
poorly understood processes which occurred dur- 



ing and soon after the birth of the star, and pos- 
sibly later due to the dynamics, it appears likely 
that different stars have envelopes with different 
chemical compositions. 

The Vela pulsar 0833-45, and possibly, but 
with much larger uncertainties, PSR 1706-44, are 
marginal candidates for enhanced cooling as their 
inferred temperatures and luminosities are lower 
than most of our models. Nevertheless, we found 
that low mass neutron stars, ~ 1.1 — 1.2 Mq , 
with a light element envelope and extensive nu- 
cleon pairing covering essentially the entire core, 
could reach the inferred values, Too — lO"''-^ K at 
an age of ~ 20,000 years (Figure 29) which is 
within the range of estimated age of the associ- 
ated supernova remnant. 

An essential component of the minimal cooling 
scenario is neutrino emission by the PBF process, 
which leads, in the presence of appropriate nu- 
cleon gaps, to more rapid cooling than possible in 
the standard cooling scenario. The low observed 
temperatures of the two pulsars. Vela and 1706-44, 
can be accomodated by the PBF process with a n 
gap of sufficient size. However, if this gap were 
vanishingly small (Schwenk & Friman 2004), then 
the temperature and thermal luminosity measure- 
ments of these two objects would be evidence for 
the presence of processes beyond the minimal cool- 
ing paradigm. 

The five older stars, PSR 0656-^14, PSR 1055- 
52, PSR 0633-1-1748 (Geminga), and RX J0720.4- 
3125, unfortunately have such large uncertainties 
on both their ages and thermal luminosities that 
their interpretation is delicate. They do not re- 
quire the occurence of enhanced neutrino emission 
and can be accomodated within the minimal cool- 
ing scenario when the lowest values for T and Loo 
are chosen. But, should the upper limits for T 
and Loo prevail, they would be good candidates 
for the occurrence of some "heating mechanism" , 
i.e., dissipative processes, which inject heat into 
the star (see, e.g., Umeda et al. (1993) & Schaab 
et al. (1999)). In the case of RX J1856.5-3754, 
which has a much more tightly constrained age 
and luminosity, the agreement with the minimal 
scenario is excellent. 

The two objects standing apart from the other 
observed neutron stars are J0205-f 6449 (in the su- 
pernova remnant 3C58) and RX J0007.0-F7302 (in 
CTA 1). Upper limits on their luminosities are 
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well below the predicted values for any of our mod- 
els. These two objects are the best candidates, to 
date, for the necessity to go beyond the minimal 
scenario. 

Finally, the four upper luminosity limits on 
the undetected objects in shell supernova rem- 
nants (G084.2-0.8, G093.3+6.9, G127.1+0.5, and 
G315.4-2.3) recently found by Kaplan et al. (2004) 
are so low that they will constitute the strongest 
evidence for enhanced neutrino emission well be- 
yond the minimal scenario if it can be demon- 
strated that they actually correspond to neutron 
stars and not quiescent black-holes. 
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Fig. 30. — Global comparison of the upper lu- 
minosity limits for sources lacking apparent ther- 
mal emissions with cooling trajectories, satisfy- 
ing the minimal cooling scenario. The area with 
the darkest gray shading contains models with 
heavy-element dominated envelopes, whereas that 
with the intermediate shading contains models 
with ligh-element envelopes. The region with the 
lightest shading contains models with intermedi- 
ate compositions. 
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A. Observational Data 



A.l. Supernova Remants G084.2-0.8, G093.3+6.9, G127.1+0.5, &; G315.4-2.3 

These four SNR's are considered to be the product of core-collapse supernovae and are hence expected 
to contain either a neutron star or a black-hole. Nevertheless, the searches of Kaplan et al. (2004) found no 
evidence of any sort for the presence of compact objects. In case a neutron star is present, these observations 
provide us with upper limits on the thermal luminosity of the star which we take from the Figure 37 of 
Kaplan et al. (2004) and report in our Table 3. 

A.2. CXO J232327.8+584842 in Gas A 

Discovered in the first light of the Chandra observatory, this object is still enigmatic but evidence points 
toward an isolated neutron star (Mereghetti et al. 2002). We take the upper limit on Lqo from Pavlov et al. 
(2000), which results from a composite model of a hot polar cap and a warm (but barely detected) surface. 
The age is from the association with the supernova SN 1680. 

A.3. PSR RX J0205+6449 in 3C58 

A first upper limit on Too had been obtained by Slane, Helfand & Murray (2002) from the non detection of 

a thermal component in the pulsar spectrum. Analysis of a deeper Chandra observation (Slane et al 2004b) 
requires the presence of a thermal component and leads to a lower upper limit on L^o reported in Table 3. 
This value is a conservative estimate since even lower values are possible when an atmosphere model is used 
in the spectral fits. The association with the historical remnant of SN 1181 gives an age of 822 years. The 
pulsar spindown age is about 5400 years (Murray et al. 2002). The distance estimate is from 21 cm (HI) 
absorption (Roberts et al. 1993). 

A.4. PSR J1124-5916 in G292. 0+1.8 

This x-ray and radio pulsar is scon as a point source in a composite SNR. The Chandra spectrum of the 
pulsar (Hughes et al. 2003) is adequately fit by a power law, with no evidence for thermal emission, and 
provides an upper limit on the thermal luminosity. Distance and kinematic age estimates are taken from 
Camilo et al. (2002). 

A.5. PSR RX J0822-4247 in Puppis A 

The ages are taken from spin-down (8000 years) and from motions of filaments in the Puppis A supernova 
remnant. The distance is estimated from 21 cm (H) absorption. Both quantities are discussed in Zavlin, 
Triimper & Pavlov (1999) and references therein. Blackbody and H atmosphere spectral fitting, from ROSAT 
and ASCA data, are also from Zavlin, Triimper & Pavlov (1999). 

A.6. PSR IE 1207-5209 in G296.9+10.0 

Estimates of the kinematic age are from Roger et al. (1988) and the distance are from Giacani et al. (1988). 
The spindown age is from Pavlov et al. (2002). Blackbody spectral fitting is from Mereghetti, Bignami & 
Caraveo (1996) from ROSAT and Zavlin, Pavlov & Triimper (1998) from ROSAT -|- ASCA. 

This is one of the very few isolated neutron stars which shows spectral lines in its spectrum (Sanwal et al. 
2002), but the phase variation of these lines (Mereghetti et al. 2002) may indicate they are of magnotespheric 
origin. Moreover, its peculiar spin-down behavior may be a sign of accretion, making the interpretation of 
this star as an isolated cooling neutron star questionable (Zavlin, Pavlov & Sanwal 2004). 

A.7. RX J0002+6246 

Spectral fitting for both H atmosphere and blackbody surfaces are from Chandra observations (Pavlov 
2002), and are consistent with ROSAT observations reported by Hailey & Craig (1995). 
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A.8. RX J0007.0+7302 in CTA 1 



An x-ray point source within a pulsar wind nebula in a composite SNR. No pulsations are (letc;c;te(l but the 
general morphology of the object makes it very similar to Vela, including a very likely association with the 
EGRET 7-ray source 3EG JOOlO+7309. The x-ray spectrum of the point source, from either XMM-Newton 
(Slane et al. 2004a) or Chandra (Halpern et al. 2004), is fitted by a power law plus a thermal component. 
The latter could be either originating from a hot polar cap or the cooler entire stellar surface. The upper 
value for the thermal luminosity we use is from Halpern et al. (2004), because the high angular resolution 
of Chandra allows a better estimate than the XMM-Newton data. Distance and kinematic age estimates are 
taken from Slane et al. (2004a). 

A.9. PSR 0833-45 (Vela) 

Spectral fitting from Chandra for both H atmosphere and blackbody surfaces are from Pavlov et al. (2001), 
and the spindown age is also quoted there. The kinematic age of the SNR is from Aschenbach, Egger & 
Triimper (1995) and the VLBI interferometric distance measurement of 250 pc is due to Dodson et al. (2003). 
This is probably the most reliable data point available and the first isolated neutron star whose radius is 
well determined because of a well known distance (Page, Shibanov & Zavlin 1996). 

A.IO. PSR 1706+44 

Blackbody spectral fitting for the Chandra data is from Gotthelf, Halpern & Dodson (2002), and spindown 
age is quoted in the same source whereas H atmosphere spectral fitting for the XMM-Newton data is from 
McGowan et al (2004). Estimates of distance are from Taylor & Cordes (1993) and Koribalski et al. (1995). 

A.ll. PSR 0538+2817 

We employ the results of spectral fitting from Chandra by Marshall & Schulz (2002) for both blackbody 
and magnetized hydrogen models. We assume a typical error of 0.1 on log^g^ of 0.5 on log^QL, since 
these authors do not report any uncertainty estimates. 

A.12. PSR J0154+61 

We use the results of an XMM-Newton observation by Gonzalez et al. (2004) for Loo- The distance and 
tsd values are taken from the same paper. 

A.13. PSR 0656+14 

We employ the results of spectral fitting from Chandra (Marshall & Schulz 2002), which result in a lower 
temperature than the ROSAT value log^^Q = 5-96l°o3^ quoted by Possenti, Mereghetti & Colpi (1996). 
Marshall & Schulz (2002) also suggest there is a hard component with a temperature of 2 x 10^ K. The 
spindown age is from Taylor, Manchester & Lyne (1993). As with other objects in this study, we employ 
the softer component's temperature as being more characteristic of the underlying surface temperature. The 
distance to this object is constrained by the VLB A parallax measurement of Brisken et al. (2003). 

A.14. PSR 0633+1748 (Geminga) 

Blackbody spectral fitting with ROSAT data is from Halpern & Wang (1997); later analyses have not 
changed these results significantly. The distance is the result of parallax measurements by Caraveo et al. 
(1996); however, Pavlov (2002) suggests these measurements may not be reliable. 

A.15. PSR 1055-52 

Pavlov (2002) quotes results from Chandra observations of two thermal components: a soft component 
with temperature 8.9 + 0.01 x 10^ K and emitting radius 13rfiooo km and a hard component with temperature 
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1.9 ± 0.1 X 10^ K and emitting radius of 0.5 ± O.ldiooo km. These tcmpcratmes arc consistent with ASCA 
results quoted in Greiveldinger et al. (1996) and ROSAT results in Ogelman (1995). We employ the soft 
component temperature as being characteristic of the average surface temperature. 

A.16. RX J1856-3754 

Absence of spectral lines in the high resolution Chandra LETGS data rules out non-magnetic, non-rotating 
heavy element atmosphere models (Burwitz et al. 2001). Blackbodies provide the best fits to the x-ray data 
but the optical data require the presence of a colder blackbody component (Pons et al. 2002). We take the 
blackbody spectral fitting from Drake et al. (2002) and Burwitz et al. (2003) for the warm component and 
fits for the colder component from Pons et al. (2002). The range of Too listed in Table 2 correspond to these 
two, cold and warm, components. The cold blackbody component gives the lower bound on the radius but 
it makes a very small contribution (~ 5%) to the luminosity. 

A.17. RX J0720.4-3125 

As for RX J1856-3754, fit of both the x-ray and optical data requires a two blackbody model. We take 
both warm and cold blackbody fits from Kaplan et al. (2003) which gives us the range of Too we report in 
Table 2. The luminosity has a significant contribution from the cold component. The distance is unknown 
and the values we report are a guess based on the observed low absorption of the x-ray spectrum. This 
distance results in a large uncertainty in Loo • -Notice also that the spectrum is known to vary on long time 
scales (de Vries et al. 2004) and contains a phase-dependent absorption line (Haberl et al. 2004). Both P 
and P are from Cropper et al. (2004). 



B. The Equations of Structure cind Evolution 

We employ the standard structure equations derived from spherically symmetric, general relativistic, 
considerations. It should be mentioned that the stellar surface in our computation is fixed by 

R = -Rstar = r{p = Pb) , (Bl) 

where = 10^" g cm^'"*. This guarantees that the EOS is temperature independent. The layers at densities 
below pb, called the envelope, are treated separately (see § 4). 

At the temperatures of interest here, neutrinos have a mean free path much larger than the radius of 
the star (Shapiro & Teukolsky (1983)) and thus leave the star once they are produced. Energy balance 
arguments (see for instance Thorne (1966)) then imply 

^ - ' 'e*(g,-gfc) , (B2) 



dr - 2Gm/c^r \dt 

where $ is the gravitational potential, I is the internal luminosity, e is the internal energy per unit volume, 
and qh are respectively the neutrino emissivity and heating rate, both per unit volume. The corresponding 
inner boundary condition for I is 

l{r = 0) = (B3) 

The time derivative of e can be written in the form 

de de dT dT 

where T is the local temperature and c„ is the specific heat per unit volume at constant volume, which for 
degenerate matter is the same as the specific heat at constant pressure cp. 
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The energy transport equation is 



dr A 47rrVl - 2Gm/c2r ' 

where A is the thermal conductivity. (Notice that within the relativistic framework, an 'isothermal' con- 
figuration is defined by e"^ ■ T = constant, instead of T = constant.) The associated boundary condition 
is 



n = n{h) , (B6) 

which relates the tcmpcratiirc T], at the outer boundary (defined more precisely further below) to the lu- 
minosity lb in this layer. The location of this outer boundary layer is chosen such that lb is equal to the 
total photon luminosity of the star, lb = l{r = R) = L. L is commonly expressed through the "effective" 
temperature Te, which is defined by 

L = AwR^-asnT^ (B7) 

where ass is the Stefan-Boltzmann constant. We emphasize that L and Tg are, modulo R, essentially 
equivalent quantities. We can thus write equation (B6) as Tb = T{,(Te); this 'Te — Tb relationship" is 
discussed further in § 4. 

We will present our results of thermal evolution by using the "effective temperature at infinity" 

r~ = Te • e*(^) (B8) 

related to the "luminosity at infinity" L°° through the "radiation radius" R^o = R • e"*^^^ by 

L°° = e2*WL(ii) = AnR°° ^ ■ assT^ ^ (B9) 

The three quantities T^,L°° and R°°, are, in principle, measurable. In particular, Roo would be the areal 
radius of the star that an observer 'at infinity' would measure with an extremely high angular resolution 
instrument (Page 1995). 
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